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ABSTRACT 

The issue of local gauge invariance in the simplicial lattice formulation of gravity is examined. We exhibit 
explicitly, both in the weak field expansion about flat space, and subsequently for arbitrarily triangulated 
background manifolds, the exact local gauge invariance of the gravitational action, which includes in general 
both cosmological constant and curvature squared terms. We show that the local invariance of the discrete 
action and the ensuing zero modes correspond precisely to the diffeomorphism invariance in the continuum, 
by carefully relating the fundamental variables in the discrete theory (the edge lengths) to the induced metric 
components in the continuum. We discuss mostly the two dimensional case, but argue that our results have 
general validity. The previous analysis is then extended to the coupling with a scalar field, and the invariance 
properties of the scalar field action under lattice diffeomorphisms are exhibited. The construction of the 
lattice conformal gauge is then described, as well as the separation of lattice metric perturbations into 
orthogonal conformal and diffeomorphism part. The local gauge invariance properties of the lattice action 
show that no Fadeev-Popov determinant is required in the gravitational measure, unless lattice perturbation 
theory is performed with a gauge-fixed action, such as the one arising in the lattice analog of the conformal 
or harmonic gauges. 
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1 Introduction 



In the quantization of gravitational interactions one expects non-perturbative effects to play an important 
role 0. One formulation available for studying such effects is Regge's simplicial lattice theory of gravity 
It is the only lattice model with a local gauge invariance and the only model known to contain 
gravitons in four dimensions [Q . One would hope that a number of fundamental issues in quantum gravity, 
such as the existence of a non-trivial ultraviolet fixed point of the renormalization group in four dimensions 
and the recovery of general relativity at large distances, could in principle be addressed in such a model. 
The presence of a local gauge invariance, which is analogous to the diffeomorphism group in the continuum, 
makes the model attractive as a regulated theory of gravity , while the existence of a phase transition in 
three Q and four dimensions ^ 0, ||, ||, |l^ (but not in two [|ll| ) suggests the existence of a (somewhat 
unusual) lattice continuum limit. The two phases of quantized gravity found in can loosely be described 
as having in one phase {G < Gc, rough, polymer-like phase) 

{g^u) = , (1.1) 

and in the other phase (G > Gc, smooth phase), 

{Qt^v) ~ crjf,^ , (1.2) 

with a small negative average curvature (anti-DeSitter space) in the vicinity of the critical point at Gc- A 
physically similar two- phase structure was later proposed also in iQ; see also the ideas found in |T^. A 
discussion of the properties of the two phases characterizing four-dimensional gravity, and of the associated 
critical exponents, can be found in For additional recent numerical results we refer the reader to pO[ |, 
while for some earlier attempts we refer to the work in [§, 0, |^ . Recently calculations have progressed to 
the point that a first calculation of the Newtonian potential from the correlation of heavy particle world 
lines, following the suggestive proposal of |l^, seems feasible ||l^. The results so far indicate that in the 
lattice quantum theory of gravity the potential between heavy spinless bodies is attractive, and has roughly 
the correct heavy mass dependence. In the same work a general scaling theory for gravitational correlations, 
valid in the vicinity of the fixed point, was put forward. We also refer the reader to p6[ |, where a more 
complete set of references to earlier work on Regge gravity can be found. For results with an alternative and 
complementary approach based on dynamical triangulations, we refer the reader to the references in [ pTf . 

In view of this recent progress it would seem desirable to further elucidate the correspondence between 
continuum and lattice theories. The weak field expansion is available to systematically develop this cor- 
respondence, and it is well known that such an expansion can be carried out in both formulations. Not 



unexpectedly, it is technically somewhat more complex in the lattice theory due to the presence of addi- 
tional vertices, as happens in ordinary lattice gauge theories. In the past most perturbative studies of lattice 
gravity have focused on the lowest order terms, and in particular the lattice graviton propagators |^. 
Recently it has been extended to include the vertex functions, and the results have been used to compute 
the one-loop amplitudes relevant for the conformal anomaly in two dimensions jl^ . 

One central issue in a regularized theory of quantum gravity is the nature of its invariance properties. 
Although some discussions of these issues have appeared before, no systematic and coherent exposition 
has been presented yet in the literature. In this paper we address the question of what exactly the local 
gauge invariance built into Regge's simplicial gravity looks like. Its existence is intimately tied in with the 
appearance of gravitons (in four dimensions) in the lattice weak field expansion about a flat background. It 
need not be emphasized here that local gauge invariance plays a central role in both the classical and quantum 
formulation of gravity, and its preservation in the lattice theory must therefore be considered of paramount 
importance. Physically, it expresses the fact that the same physical geometry can be described by equivalent 
metrics. Classically, it leads for example to the invariance of the infinitesimal line element and the Bianchi 
identities for the curvature. In the quantum theory it is known to give rise to the Slavnov- Taylor identities 
for the gravitational Green's functions. One would therefore expect that local gauge transformations should 
play a central role in the lattice theory as well. This aspect will be therefore the focus of the first part of the 
paper, where the analog of local gauge transformations on the lattice will be constructed. The requirement 
of gauge invariance will have implications for both the gravitational measure and the coupling to a scalar 
field, and we will present in this paper a detailed analysis of its consequences. The second part of this paper 
will be devoted to a number of relevant applications. 

The plan of the paper is as follows. In Section 2, we introduce our notation, describe the choice of 
lattice structure and the relevant degrees of freedom in the lattice theory, the squared edge lengths. We 
discuss the discrete actions for the gravitational degrees of freedom, and the relationship to their well- 
known continuum counterparts. In Section 3 we move on to the lattice weak field expansion, and discuss 
in detail the two-dimensional case (with cosmological and curvature squared terms). We exhibit explicitly 
the gauge zero modes and their corresponding eigenvectors, which are shown to correspond precisely to 
local gauge transformations in the continuum. We then compute explicitly and analytically the zero modes 
for fluctuations about a non-flat background (the tetrahedral, octahedral and icosahedral tessellations of 
the two-sphere), and show that the counting of the zero modes is consistent with the expectation from the 
continuum theory We then give further arguments supporting the identiflcation of the zero modes with 
the diffeoniorphisms in the continuum, which we argue is valid in any dimension. In Section 4 we extend 
the previous analysis to arbitrary curved backgrounds and show explicitly the persistence of a local gauge 



invariancc for the area, curvature and curvature squared terms. In section 5 we introduce a scalar field 
coupled invariantly to the gravitational degrees of freedom. We again exhibit its invariance properties under 
local gauge variations of the squared edge lengths, at least for sufficiently smooth scalar field configurations, 
by working out the concrete case of background lattices which are close to either equilateral or square. We 
then discuss the more general case of arbitrary background lattices, and the construction of the energy- 
momentum tensor for the scalar field. Section 6 discusses the implications of the preceding results for the 
lattice gravitational measure, and we give arguments that the lattice measure is essentially unique, up to 
local volume factors. We will argue therefore that the lattice measure is essentially no less unique than 
the original continuum (DeWitt) measure. In Section 7 we consider the possibility of introducing a gauge 
fixing term in the lattice action, in order to remove the gauge zero modes of the gravitational action and 
subsequently perform perturbative calculations, and in close analogy with the procedure followed in the usual 
continuum perturbation theory. As an example, we discuss the explicit construction of the lattice conformal 
gauge, starting from an arbitrary configuration of squared edge lengths. Finally, Section 8 contains some 
concluding remarks. 



In this section we will briefly review the construction of the action describing the gravitational field on 
the lattice, and define the necessary notation used later in the paper. In concrete examples we will often 
refer, because of its simplicity, to the two-dimensional case, where a number of results can be derived easily 
and transparently. But in a number of instances here, and throughout the paper, important aspects of 
the discussion and of the conclusions will be quite general, and not restricted to specific aspects of the 
two-dimensional case. 

In simplicial gravity the elementary building blocks for rf-dimensional space-time are simplices a'^ of 
dimension d. A 0-simplex is a point, a 1-simplex is an edge and a 2-simplex is a triangle. A rf-simplex is 
a d-dimensional object with + 1 vertices and d{d + l)/2 edges connecting them. Each simplex in turn 



closed simplicial two-manifold consisting of N^y vertices, iVi edges and N2 triangles, joined in such a way 
that each point has a neighborhood homeomorphic to the interior of a two-dimensional sphere. A simplicial 
geometry is then specified by the assignment of squared edge lengths lf,i = l... Ni, and a flat Riemannian 
metric can be assigned to the interior regions of the simplices in a way that is consistent with the edge length 
values. Further restrictions arise from the fact that the triangle inequalities (and their higher dimensional 
analogs in d dimensions) have to be satisfied. 



2 The Discretized Theory 




Thus in two dimensions we shall consider here a fixed 



The correspondence between squared edge lengths and an assigned continuum metric field can be made 
more precise, with the identification 

r(b) I . 



where lab is the length of the edge connecting neighboring points a and b. For a given set of edge lengths, the 
metric g^u{x) has initially support on the edges only. ^ For a metric that is constant inside each simplex. 
If = g^ul^l^, where i labels the edge from a to 6 and the 's are the components of the edge lengths. 

2.1 Lattice Structure 

In two dimensions quantum gravity can be defined on a two-dimensional surface consisting of a network 
of flat triangles. The underlying lattice may be constructed in a number of ways. Points may be distributed 
randomly on the surface and then joined to form triangles according to some algorithm. In such lattices the 
coordination number at each vertex can be kept fixed (quenched random lattice) , or allowed to vary (annealed 
random lattice), by considering it as an additional, dynamical variable of the model. An alternative procedure 
is to start with a regular lattice, like a regular tessellation of the two sphere, or a lattice of squares divided 
into triangles by drawing in parallel sets of diagonals, and then allow the edge lengths to vary, which will give 
rise to curvature localized on the vertices. It should be emphasized that for arbitrary assignments of edge 
lengths, consistent with the imposition of the triangle inequalities constraints, such a lattice is in general far 
from regular, and resembles more a random lattice. 

The incidence matrix, which provides the information on which edges are adjacent, and fixes therefore 
the local coordination number g^, describes the topology of the manifold. It can be chosen to correspond to a 
fixed regular or to a fixed random lattice. But one word of caution should be spent here on the terminology. 
Since the edge lengths are dynamical variables, the lattice is in fact random in either case: contrary to a 
fixed regular lattice (such as the square or triangular one in two dimensions), there are a priori no preferred 
directions even for a lattice with fixed coordination number, as neighboring points can have any relative 
orientation as long as they are consistent with the triangle inequalities and their higher dimensional analogs. 
Universality arguments would then suggest that the choice of local coordination number should not affect 
the large distance limit of the model, and a number of explicit calculations on random lattices have shown 
to some extent that this is indeed the case pO[ . 

^The above identification parallels an analogous correspondence used sometimes in ordinary lattice gauge theories, 
where the SU{n) matrix- valued lattice field Anfi has support only on the links of a hypercubic lattice, Unfi = e*°'^"'' = 

Pexp ^ia J'"^'' dxf^Afj_(x)^ . This definition is a convenient starting point for performing perturbation theory and defining the 

lattice Feynman rules. For the same construction in pcrturbative simplicial gravity see [|L8|. 




(2.1) 



In the following we will often narrow down the discussion and be even more specific, and usually think of 
the "regular" lattice as consisting of a network of triangles with a fixed coordination number of six, qi =6, 
although many of the results in this work are quite general and do not rely on the specific choice of local 
coordination numbers. 

Quenched random lattices, where the local coordination number qi (which is the number of edges meeting 
at i) is random but fixed, were considered in |2^, [l^, For such Poissonian random lattices, the average 
coordination number is also q — 6 in two dimensions, independent of the topology. This follows from the 
expression for the Euler characteristic x — Nq ~ Ni + N2 with 2Ni — 3N2 — J2i 1i two dimensions, 
which gives for large A'o 

q^ \im ^ = 6 , (2.2) 

irrespective of the value of x, as well as iVi = SNq and N2 — 2iVo. In general on such random lattices one 
does not have, strictly speaking, translational or rotational invariance for a fixed assignment of edge lengths. 
The latter only hold on the average. Explicit calculations confirm that this is indeed the case, at least in 
two dimensions |^ |5|. 

When the edge lengths are allowed to fluctuate one would expect the situation to be different, since 
now locally there are no preferred directions any more, as the lattice structure fluctuates from edge length 
configuration to edge length configuration. Lastly, one can allow the local coordination number to change 
(annealed random lattice) by re-linking neighboring vertices, although there is no unique algorithm to do 
so which preserves the geometry. In this case the coordination number fluctuation Sqi — qi — 6 becomes 
an additional dynamical variable, and is indeed the only dynamical variable in the so-called dynamical 
triangulations. Randomness can be shown to be a relevant perturbation in two dimensions, changing the 
universality class already for flat surfaces. We shall not consider dynamical random lattices here any further, 
as we are interested in discretizations for gravity coupled to matter which maintain the crucial property of 
reducing to the ordinary, known flat space field theories in the limit of zero local curvatures. A review of the 
properties of random lattices and their relation to matrix models in two dimensions can be found in [ p6[ . 

2.2 Degrees of Freedom 



The elementary degrees of freedom on the lattice are the edge lengths, with the correspondence between 
continuum and lattice degrees of freedom given locally by 

{5,.(^)L,^-{^n.=i...Ar, > (2-3) 

where the index i ranges over all N edges in the lattice. In general the dynamical lattice will give rise to 

1 

some average lattice spacing ao = , which in turn will naturally supply the ultraviolet cutoff that is 



needed to define the quantum theory. An important difference with ordinary lattice field theories lies in the 
fact that the momentum cutoff A — i/Zo is not determined a priori, but follows instead from the dynamics 
(i.e. from the lattice action and lattice measure). The dynamical cutoff turns out to be determined mostly 
by the cosmological constant term and the measure factor 

Furthermore for finite volumes the lattice theory will have a finite number of degrees of freedom N , and 
will therefore inherit an infrared cutoff of the order of l/L, where L is the physical linear extent of the 
lattice. 

In the discrete case all the metric information on the piecewise linear space is contained in the values 
of the edge lengths. As already emphasized by Regge, and in accordance with the usual view of lattice 
discretization of continuum field theories, the discrete manifold L is thought of as an approximation to some 
continuum manifold S (this is illustrated in Figure 1). In the limit as the average lattice spacing ao is sent 
to zero, the original continuum theory is recovered. In four dimensions it has been rigorously proven, for the 
Einstein-Regge action, that if a piecewise flat space approximates a smooth space in a suitable sense, then 
the corresponding curvatures are close in the sense of measures ; see also the results of . In general 
the expectation is that the lattice and continuum theory will differ by higher order corrections, with the two 
actions related to each other by 

- /c(<?M-) + ao'5/ + ao'5'/+--- • (2.4) 

All corrections can in principle systematically be evaluated by the standard procedure of replacing the finite 
differences which appear in the lattice action by derivatives, for example according to the formula 

g(n + ao)-g(n-ao) // n 1 i,i , n , ^/ 4 (5)/ /n r^ 
= g{n)^-ang (n) + O(ao g^^^(n)) . (2.5) 

It should be noted that higher order corrections are expected to involve higher derivatives of the metric. 
The above expansion procedure can be thought of being equivalent to introducing a continuum metric on 
the piecewise linear manifold, and expand in the difference between the continuum and the piecewise linear 
metric. 



Fig. 1. Piecewise linear space L as an approximation to a smooth d-dimensional enveloping surface S. 

This interpretation is analogous to the situation in ordinary lattice gauge theories, where the lattice 
gauge fields ?7„^ are defined on the links only; the continuum fields Afj_{x) can then be reconstructed by 
some suitable interpolation to the interior regions of the lattice. It is of course possible to endow the 
piecewise linear space with a continuum metric g^i,^{x) which is defined everywhere, including the interiors of 
the simplices. In this case a continuum curvature R^vpaix) can be defined as well, but since the interior of the 
simplices is flat, the curvature acquires delta- function singularities on the hinges where the discrete curvature 
resides. While such a description can be useful in certain circumstances, it has also some drawbacks, which 
have led to considerable confusion in some of the literature. The obvious ones are that the resulting model 
is no longer an ultraviolet regulator for the continuum theory, as space-time has become continuous again. 
Furthermore the fields are singular, due to the delta-function singularities in the curvature, and the number 
of degrees of freedom is no longer finite due to the re-introduction of a continuum metric. In this formalism 
new divergences appear, which have to be regulated by some ad-hoc procedure such as the smoothing out 
of conical singularities, and lead to difficulties in defining higher order invariant operators such as the ones 
containing curvature squared terms 0] . 

This point of view, while certainly legitimate in discussing some classical aspects of the theory, is therefore 
in our opinion not useful in describing a regulated theory of quantum gravity. It leads instead to a string of 
paradoxical results when lattice and continuum language are mixed together, and can be especially misleading 
when discussing such subtle issues as the gravitational functional integration measure. 

It should be emphasized here that in the following we shall restrict our attention almost exclusively to 
the lattice theory, which is defined in terms of its lattice degrees of freedom only. Since it is our purpose to 
describe an ultraviolet regulated theory of quantum gravity, we shall follow the usual procedure followed in 
discussing lattice field theories, and discuss the model exclusively in terms of its primary, lattice degrees of 
freedom: the squared edge lengths. As such, the theory will not require any additional ad-hoc regulators. 
Below we shall discuss further at length a number of issues related to the precise correspondence between 
the lattice degrees of freedom and the continuum ones, the local gauge invariance of the lattice action (which 
gives rise in the quantum theory the lattice analogs of the Taylor-Slavnov identities) and the need for (or 
lack of) gauge fixing. 

2.3 Curvature and Discretized Action 

The construction of the lattice action starts from the definition of the elementary building blocks for 
spacetime, the n-dimensional simplices. Consider an n-dimensional simplex with vertices 1, 2, 3, ... n-l-1 



and square edge lengths If^ =■ %\i ••• • Its vertices are specified by a set of vectors eb = 0, ei, ... e„ in flat 
Euclidean space. The matrix 

ga = ei ■ Cj , (2.6) 
with 1 < i,i < n, is positive definite. In terms of the edge lengths = |e, — e*, ! (see Figure 2) it is given by 

9ijil') = \ + ■ (2.7) 

The volume of a general n-simplex is then given by an n-dimensional generalization of the well known formula 
for the volume of a tetrahedron, 

Vn{e) = ^^det5i,(/2) . (2.8) 

Conversely, in order to obtain a simplex for an arbitrary assignment of edge lengths, the generalization to 
higher dimensions of the triangle inequalities require that V^'^ (Z^) > 0, with n = l . . .d and i = 1 . . . Nn be 
satisfied for every edge, triangle, tetrahedron etc. in the lattice. This can be stated equivalently by requiring 

detgij{l^) > (2.9) 

for every sub-determinant of the highest dimension detg^ . In d dimensions the matrix gij has d(d + l)/2 
components, just as there are d{d + l)/2 components for the metric g^v{x) per space-time point in the 
continuum. 




Fig. 2. Assignments of edge lengths for a four- dimensional simplex. In this paper we shall often refer to the 
two-dimensional case. In two dimensions one has simply 

n..(i'^\ — ( 'oi ^('oi + 'o2 ^ '12) \ O ^(\^ 

V 2 v'Ol "T '■02 '■12) '02 / 



and therefore 

det gij{P) — I [2(^01^02 + ^02^12 + ^12^01) ~ ^01 ^ ^02 ~ ^12] (2-11) 

and 



yciet5~(P) = 2AT{f) , (2.12) 
where At{P) is the area of the given triangle (see Figure 3). 




Fig. 3. Assignments of edge lengths and natural coordinates for a triangle. 

In simplicial gravity the curvature is concentrated on the hinges, which are subspaces of dimensions d — 2, 
and is entirely determined from the assignment of the edge lengths. In two dimensions the hinges correspond 
to the vertices and 6h, the deficit angle at a hinge, is defined by 

6h = 27:- J2 ^* ' (2.13) 

triangles t 
meeting at h 

where dt is the dihedral angle associated with the triangle t at the vertex h (see Figure 4). In d dimensions 
several d-simplices meet on a (d — 2)-dimensional hinge, and the deficit angle is defined by 

Sh{l^) = 2n- g 0^{l^) , (2.14) 

meeting on h 

where 9d is the dihedral angle in d dimensions. The sine of the dihedral angle can be computed from the 
well known formula 

sin9d{l ) - — -J— , (2.15 

d-1 Vd-iV^_i 

where Vd-2 is the volume of the hinge, Vd is the volume of the d-simplex, and Vd-i, V^_i the volumes of 
the two {d — l)-dimcnsional faces that meet on the hinge. A general derivation of these formulae can be 
found in with some additional results in 0. Since the sine does not uniquely determine the angle, it 
can be useful to obtain an expression for the cosine of the dihedral angle, which can be found in 0. In two 
dimensions the dihedral angle is given by 



Fig. 4- In two dimensions the computation of the deficit angle 5o at the vertex involves the values for the edge 
lengths associated with the shaded triangles. 

It is useful to introduce a dual lattice following, for example, the Dirichlet-Voronoi cell construction, which 
consists in introducing the perpendicular bisectors of the edges in each triangle and joining the resulting 
vertices. This provides for a natural subdivision of the original lattice in a set of non-overlapping exhaustive 
cells, and has a natural generalization to higher dimensions. It is easy to see that the vertices of the original 
lattice then reside on circumscribed circles, centered on the vertices of the dual lattice. For the vertex the 
dihedral dual volume contribution, shown in Figure 5, is given by 

(2.17) 

It is clear that the above subdivision is not unique. Alternatively, one can introduce a baricenter for each 
triangle, defined as the point equidistant from all three vertices, and again join the resulting vertices. The 
vertices of the original lattice then reside on inscribed circles, centered on the vertices of the dual lattice. 
The baricentric dihedral volume is simply given by 

Aail^) = A/i . (2.18) 

In general, if the original lattice has local coordination number at the site i, then the dual cell centered 
on i will have qi faces. A fairly complete set of formulae for dual volumes relevant for lattice gravity and 
their derivation can be found in 0. In the following we shall refer to the Voronoi cell construction as the 
"dual subdivision" , while we will call the baricentric cell construction the "baricentric subdivision" . 



^12(^01+^02) (^01 ^02)^ 



2 




Fig. 5. Dual area Ad associated with vertex 0, and the corresponding dihedral angle Od- 

Two-dimensional Einstein gravity is trivial because the Einstein action is constant and the Ricci tensor 
vanishes identically. When a cosmological constant term and a curvature-squared term are included in the 
action, 

(fxyr^\x-kR + aR^] , (2.19) 



the classical solutions have constant curvature with R = ±^/xJa (there being no real solutions for A < 0). 
Thus the theory with the Einstein action and a cosmological constant is metrically trivial, having neither 
dynamical degrees of freedom nor field equations. On the other hand the functional measure can lead to a 
non-trivial effective action. However, for a system with fixed topology, the only non-classical aspects of 1 + 1 
dimensional gravity are fluctuations in the local volumes ^y g{x). 

The Einstein action for a two-dimensional simplicial lattice is given by 

/ d^x^R -^2 J2 ■ (2.20) 

hinges h 

According to the Gauss-Bonnet theorem the Einstein action in two-dimensions is equal to 47r times the Euler 
characteristic of the surface. The same result is true on the lattice, with ~ 27rx, where x is the Euler 

characteristic. It is a constant provided we consider, as we shall do below, surfaces with a fixed topology. 
A cosmological constant term can be included in the action in the form 

A j d?x^ \ J2 ^* ' (2.21) 

triangles t 

where At is the area of triangle t. Equivalently one may subdivide the triangles into areas associated with 
each hinge Ah and use the expression 

X Ah . (2.22) 

hinges h 



For the baricentric subdivision one has simply 



Ah - 



triangles t 
nccting at 



(2.23) 



Ah can also be taken to be the area of the cell surrounding h in the dual lattice (see Figure 6) , with 



Ah 



E 



A, 



(2.24) 



trianglcB t 
meeting at h 



with the dual area contribution for each triangle given in Eq. ( p. 17 ) 




Fig. 6. Original simplicial lattice (continuous lines) and dual lattice (dotted lines) in two dimensions. The shaded 
region corresponds to the dual area associated with vertex 0. 

In two dimensions the Weyl tensor vanishes identically, while the other curvature-squared terms are all 
proportional to each other, 

RnypaR'^"'"^ — \ RpiyR^^ = R? ■ (2.25) 

One therefore needs only one term quadratic in the curvature for the lattice action. Using the requirements 
that it be a sum over hinges (the only places where the curvature is non-zero), that it be quadratic in the 
deficit angle, and that it have the correct dimension (length) one is led to the unique expression 

j<fx^gR^ ^aY,^^ . (2.26) 

It can be shown that this formula is exact for all regular tessellations of the two-sphere, in the sense that the 
discrete lattice expression does not depend on the how fine the tessellation is, once the area of the surface is 
kept fixed Q. 



The lattice action corresponding to pure gravity is then 



, L Ah 



which can be written equivalently as 



(2.27) 



(2.28) 



with the two-dimensional volume element Vh — Ah, and the local curvature given by Rh = 2 Sh/Ah- In the 
limit of small fluctuations around a smooth background, corresponds to the continuum action 



m 



\-kR + aR' 



(2.29) 



For a manifold of fixed topology the term proportional to k can be dropped, since Sh — 27rx, where x is 
the Euler characteristic. The curvature-squared leads to non-trivial interactions in two dimensions, although 
the resulting theory is not unitary. In the next section we shall discuss properties of the above action in the 
weak field expansion about flat space, and later about an arbitrary lattice manifold. 

Arguments based on perturbation theory about two dimensions (where the gravitational coupling is 
dimensionless and the Einstein theory becomes renormalizable) suggest that there should be no non-trivial 
ultraviolet fixed point of the rcnormalization group in two dimensions. Explicit calculations in the lattice 
theory have shown conclusively that this is indeed the case in the absence of matter |o[|5, ||, || . The 
equations of motion for pure gravity in two dimensions then follow from the variation 



SI[g] = ijd'x^ [A - g'^'Jg^, - 



and read 



a 2 A 



or, in contracted form, R^ — ^. For an arbitrary gauge variation of the metric, 

^9t^,y{x) = -gf,\{x)d^x^{x) - g\,,{x)d^x^{x) - dxg^^{x)x^{x) , 



(2.30) 



(2.31) 



(2.32) 



one obtains after an integration by parts, and using the fact that the gauge function x^ is arbitrary (and 
that {gn.,. = 0), 

(i^V)^, = . (2.33) 

This is the two-dimensional analog of the (contracted) Bianchi identity. Since the squared edge lengths are 
the primary degrees of freedom, the corresponding lattice field equations of motion are obtained, in any 
dimension, from 

-- . (2.34) 



di 



Already in the two-dimensional case they are rather unwieldy when written out explicitly, and will not be 
recorded here. 

A candidate for the discrete analog of the two-dimensional Bianchi identity is simply 

Y^S,{f^+6l^) - = , (2.35) 

h(i) h{i) 

where the sum includes the four hinges h belonging to the two triangles bordering the edge i, and Slf 
represents a variation of the edges meeting at the vertex h. By considering gauge variations of the edge 
lengths in higher dimensions, the corresponding exact lattice Bianchi identities can in principle be written 
down. Some further discussion of the Bianchi identities in higher dimensions can be found in the second 
reference in [Q. 

3 Lattice Weak Field Expansion and Zero Modes 



One of the simplest problems which can be studied analytically in the continuum as well as on the lattice 
is the analysis of small fluctuations about some classical background solution. In the continuum, the weak 
field expansion is often performed by expanding the metric and the action about flat Euclidean space 

9^ll'{x) = S^^ + Kh^^{x) . (3.1) 

In four dimensions k = vMttG, which shows that the weak field expansion there corresponds to an expansion 
in powers of G. In two dimensions this is no longer the case and the relation between k and G is lost; instead 
one should regard k as a dimensionless expansion parameter which is eventually set to one, k = 1, at the end 
of the calculation. The procedure will be sensible as long as wildly fluctuating geometries are not important 
in two dimensions (on the lattice or in the continuum). The influence of the latter conflgurations can only 
be studied by numerical simulations of the full path integral |jri|, ^ . 

In the lattice case the weak field calculations can be carried out in three |^ and four |Q] dimensional flat 
background space with the Regge-Einstein action. One finds that the Regge gravity propagator indeed agrees 
exactly with the continuum result |Q in the weak-field limit. As a result, the existence of gravitational 
waves and gravitons in the discrete lattice theory has been established (indeed it is the only lattice theory 
of gravity for which such a result has been obtained ^ ) . 

The weak field expansion about flat space is relevant for the continuum limit of the lattice quantum theory. 
Consider a simplicial lattice approximation to a given continuum manifold. For an arbitrary continuum 

*A discretization of the edge lengths, and therefore of the curvatures, as advocated in some models for lattice gravity, 
can be considered, the dynamical triangulations being one specific example. This procedure leads obviously to a loss of the 
graviton excitation, at least in the weak field expansion. In ordinary non-abelian lattice gauge theories, models based on discrete 
subgroups of SU (N) have an artificial freezing transition at finite c oup ling and no lattice continuum limit, and do not seem to 
represent a useful discretization of the original continuum theory [p5| . 



manifold, one can envision a triangulation which is successively refined by making the simplices and the 
corresponding edge lengths smaller and smaller. As the average lattice spacing is reduced, the curvature 
on the scale of the lattice spacing becomes eventually sufficiently small that the simplicial manifold can be 
regarded as being locally close to flat. In this limit the curvature is small on the scale of the local volume, 
and in two dimensions one has 



\curvature\h 



5h 



< {volume),^^ = \- or < 1 . (3.2) 



In such regions, which become larger and larger in size as the lattice spacing is reduced, one can meaningfully 
apply the weak field expansion about flat space, which becomes only an approximation when it is truncated 
to any finite order. 

In the following we shall consider in detail only the two-dimensional case, although similar calculations 
can in principle be performed in higher dimensions, with considerable more algebraic effort. In the pure 
gravity case the Einstein-Regge action is a topological invariant in two dimensions, and one has to consider 
the next non-trivial invariant contribution to the action. We shall therefore consider a two-dimensional 



lattice with the higher derivative action of Eq. (2.27) and A = 0, 

m = ia Y: ^ • (3.3) 

hinges h 

The weak field expansion for such a term has largely been done in |pj]| , and we will first recall here the main 
results. Since fiat space is a classical solution for such an i?^— type action, one can take as a background 
space a network of unit squares divided into triangles by drawing in parallel sets of diagonals (see Figure 7). 
This is one of an infinite number of possible choices for the background lattice, and a rather convenient one. 
Physical results should in the end be insensitive to the choice of the background lattice used as a starting 
point for the weak field expansion. Opposite edges of the network are supposed to be identified so that the 
lattice acquires the topology of a torus. (In the following we will be concerned with local properties of the 
action, and the detailed nature of the boundary conditions will play only a marginal role). 



(0,1) 




Fig. 7. Notation for the weak-field expansion about the rigid square lattice. 

It is also convenient to use the binary notation for vertices described in references [Q. As discussed 
in the previous section, the edge lengths on the lattice correspond to the metric degrees of freedom in the 
continuum. The edge lengths are thus allowed to fluctuate around their flat space values, 

k = 1° (1 + e,) , (3.4) 

with If^l^^l and ?g = V2 for our choice of lattice. The second variation of the action is then expressed 
as a quadratic form in the e's, 

S'^I ^ Aa^ei M,j . (3.5) 

The properties of Afy are best studied by going to momentum space. One assumes that the fluctuation 
at the point i, j steps in one coordinate direction and k steps in the other coordinate direction from the 
origin, is related to the corresponding Ci at the origin by 

where oji = e~'^^' and ki is the momentum in the direction i. The matrix M then reduces to a 3 x 3 matrix 
with components given by 

{M^)ii = 2 + ivi — 2uj2 - 2iLjiUj2 + ojiiol + c.c. 

(Mt^)l2 = 2 — UJi — u)2 — ^1^2 — ijJlOj2 — ^ + + ijJlijj2 + 20^1^2 

{M^)lZ = 2{—l + 2lUi — LUi + 1^2 — ^2 — l^l^2 + '21^1^2 + ^2 ^ ^1^^2 ^ 1^1(^2) 

{Muj)33 = 4(2 — 2wi — 2w2 + 0J1UJ2 + OJ1UJ2 + c.c.) 

(3.7) 

with the other components easily obtained by symmetry. For small momenta M^^ takes the form 

/ P2{ki+k2)^ kik2{k^+k2f -2kikl{ki + k2) \ 

kik2{ki + k2f kl{ki+k2f -2klk2{ki + k2) +0{k'-) . (3.8) 
\ ~2kikl{ki + k2) -2kjk2{ki + k2) Aklkl j 

The change of variables 

e'l = ei 4 = £2 4 = 2(ei+e2) + e3 • (3.9) 
leads for small momenta to the matrix M'^ given by 

M'^ = I klkl kf -2fc3fc2 I + 0{k^) . (3.10) 



-2kit^ -2fcffc2 ikfk 



2 



This expression is identical to what one obtains from the corresponding weak-field limit in the continuum 
theory. To see this, define as usual the small fluctuation fleld /i^^ about flat space by setting 

(3.11) 



In two dimensions one has 



and also 



R = /111.22 + /l22,ll - 2/ll2,12 + , (3.12) 



V5 = l + i(ftii + /i22)+0(/i') , (3.13) 



which gives 



VgR^ - (/iii,22 + Vii - 2/112,12)' + o(/i3) . (3.14) 

In momentum space, each derivative produces a factor of fc^, and so one obtains 

y/gR^ = h^yV^y^pahp^ , (3.15) 

where V^u.pa coincides with M' above (when the metric components are re-labeled according to 11 1, 22 —> 
2, 12 3). 

One might wonder what the origin of the change of variables in Eq. (|]^) is. Given the three edges in 
Figure 7, one can write for the metric at the origin 



The apparent contradiction with the earlier expression for gij given in Eq. (2.10) arises from the different 
choice of coordinates in the triangles (compare Figure 3 with Figure 7). Inserting U = (1 + 6^), with = 1 
for the body principals (i = 1, 2) and = V2 for the diagonal (i = 3), one obtains 



li = (1 + ei)' = l + hii 

ll = (1+62)2 = l + h22 
\ll = (1+63)2 = l + i(/lll+/l22) + /jl2 



(3.17) 



which can be inverted to give 



ei = lh^i-lhl,+0{hl^) 

62 = i/l22- i/l22+0(/l^2) 

£3 = i(/lll+/l22 + 2/742) -52(^11+ ^22 +2/li2)2 + 0(/l3) 



(3.18) 



and which was then used in Eq. ( |3.9D . Thus the matrix was brought into the continuum form after 
performing a suitable local rotation from the local edge lengths to the local metric components. 

The weak field expansion for the purely gravitational part can be carried out to higher order, and the 



Feynman rules for the vertices of order h^, ft,*, ... in the i?^-action of Eq. (2.27) can be derived. Since their 
expressions are rather complicated, they will not be recorded here. 

3.1 Lattice Diffeomorphisms 



It is easy to determine the eigenvalues and eigenvectors of the matrix M^^ of Eq. ( p.7[ ) . The eigenvalues 
of the matrix M^^ are given by 

Ai = 
A2 = 

A3 = 24 - 9(cJi + CJi + CJ2 + tJ2) + 4(^1^2 + <I'lt^2) 
+UJ1UJ2 + Ct;^CJ2 + 'I'l'I'l + <^Jct'2 

(3.19) 

and there are thus two exact zero modes in the weak field limit. It should be emphasized that the exact 
zero modes appear for arbitrary uji, and not just for small momenta. We shall see later that their presence 
reflects an exact local continuous invariance of the gravitational action. 

If one were interested in doing lattice perturbation theory, one would have to add a lattice gauge fixing 
term to remove the zero modes, such as the lattice analog of the term 

^ {d,y/^) g^'^y , (3.20) 

and add the necessary Fadeev-Popov non-local ghost determinant. A similar term would have to be included 
as well if one were to pick the lattice analog of the conformal gauge |3^, to which we shall return later. 
If one is not doing perturbation theory, then of course the contribution of the zero modes will cancel out 
between the numerator and denominator in the Feynman path integral representation for operator averages, 
and such a term should not be included, as in ordinary lattice formulations of Yang-Mills gauge theories. 
The eigenvectors corresponding to the two zero modes can be written as 

lit]] = i '^o"' l'., ](lfj)) ' (3.21) 

esiuj) J V 1(1-^1^2) i(l-c.ic.2) / ^ ' ' 

where Xi('-^) X2('-^) f^i'e arbitrary. One might worry that the above result is restricted to two dimensions. 

This is not the case. Completely analogous zero modes are found for the Regge action in three |0| and 



four [Q dimensions, leading to expressions rather similar to Eq. (3.21), with as expected d zero modes in 
d dimensions. As we shall see, this is not a coincidence. We give here for comparison the corresponding 
expression in three dimensions obtained from the weak field expansion of the Regge action, 

\ / l-LUl \ 



1 — U!l 




i(l - UJ1UJ2) 
i(l - UJ1UJ4) 





I-UJ2 


i(l-u;iW2) 

i(l-CJ2W4) i(l-a;2C^4) 




(3.22) 



again with xi(a;), X2{^) a-nd Xsi^) arbitrary gauge functions (in the binary notation for the edges, the 
indices 1,2,4 correspond to the body principals, the indices 3,5,6 to the face diagonals, and 7 to the body 
diagonal). 



X (n + ^j) 



X (n + Hi + lij) 




Fig. 8. Edge length gauge deformations ei{n), and corresponding gauge transformation vector field x(^)) defined on 
the sites. 



It is useful to look at the above relations, and in particular Eq. (3.21, in real space. The replacement 
Q^k, _^ e^'^ f{x) = f{x + a) yields 

ei(") = Xi{n) - Xi{n + M 
e2(n-) = X2{n) - X2{n + fi2) 

^3{n) = 5X1(^1) + 5X2(n-) - 5Xi(" + Ai +A2) - 5X2(^1 + Ai +A2) ■ 

(3.23) 

For our notation, we refer to the drawing in Figure 8. Note how the arbitrary gauge variations act on the 
two ends of an edge. This is true in any dimension, where elementary local gauge transformations are always 
defined on the vertices of the lattice. The above gauge transformation law is in fact remarkably simple for 
those edges that lie in the direction of the chosen coordinates, namely 

= X'^ - X'j : (3.24) 



where i and j labels the end-points of the edge, and we have rescaled the arbitrary functions x by Ifj so that 
the quantities x' i now have dimensions of length squared. 

It is easy to see that the above equations solve the constraint £i{n) = at the vertex n where Xn 7^ 0, 

6 

X! ^»(") Xiin) + + ^X2{n) + X2{n) - Xi(") - 5Xi(") - ^X2{n) - X2{n) 

= . 

(3.25) 

It can be written equivalently in terms of variations of the squared edge lengths meeting at the vertex n, 
labeled clockwise around the vertex n starting with the edge in the positive 1 direction, 

a {Ki + + + + b {Sll^ + Sll,) = , (3.26) 

with a and b arbitrary constants at this point. In other words, gauge variations of the squared edge lengths 
are recognized as special variations, where all edges meeting at a point (in two dimensions) are considered, 
and which either have the explicit form given in Eq. ( ^.23D for the weak field case, or equivalently (and more 



generally) satisfy a set of defining constraints such as the one in Eq. (3.26). 

Incidentally, it should be noted here that conformal transformations, which in the continuum take the 
form 5gfj_^{x) = Qiivix) Sip{x), have a natural lattice analog. They contract (or expand) locally all the edges 
meeting on a given vertex n by the same amount, 

Sl^n) = l^{n) Sipin) , (3.27) 

and therefore do change locally the curvature at n. When constructed in this way, they can be considered 



orthogonal to the x transformations of Eqs. (3.21) and (3.23). 

The cosmological constant term can also be shown to be invariant under the same set of continuous local 
transformations, since, using the same notation for the expansion about the square lattice, one obtains 

= ^[l + i(ei(n)+e2(n)) + 0(e^)] , (3.28) 

h n 

and at the vertex P„ where Xn 7^ one has again 

4 

= Xi{n) + X2{n) - Xi(n) - X2{n) = , (3.29) 

i=l 

where the sum is over the four edges pointing in the four principal directions. Written in terms of the 
variations of the squared edge lengths, one has 

+ + = . (3.30) 



In general one has in two dimensions six edges meeting at a vertex, and therefore four aUowed constraints 
on the gauge edge length variations. In conclusion we have exhibited an exact local gauge invariance of the 
gravitational action in the weak field limit. Later on we shall show that it corresponds precisely to the lattice 
analog of the diffeomorphisms. 

It is important to notice that the appearance of zero modes in the weak field expansion is not specific 
to the expansion about flat space. One can look at the same procedure for variations about spaces which 



are classical solutions for the gravitational action with a cosmological constant term as in Eq. (2.27), such 



as the regular tessellations of the two-sphere [Q. In the following we will consider edge length fluctuations 
about the regular tetrahedron (with 6 edges), octahedron (12 edges), and icosahedron (30 edges). 





Fig. 9. Tetrahedral tessellation of the two-sphere, with arbitrary edge length assignments. 

After expanding about the equilateral configuration, the action at the stationary point reduces to 

/ = A 87rVaA + a Stt/VoA = 167r\/^ , (3.31) 

in fact independently of the tessellation considered. Vanishing of the linear terms in the small fluctuation 
expansion gives for the average edge length 

lo = [cTr\Aa/X)Y^^ , (3.32) 

with c = 16/3,4/3,16/75 for the tetrahedron, octahedron and icosahedron, respectively. For fluctuations 
about the classical solution for a tetrahedral tessellation of S'^ (see Figure 9) the small edge length fluctuation 
matrix gives rise to the following coefficients 

e?2 16 VoA (54 - gVStt + 57r2)/8l7r 

£12 613 16\/aA7r/9 

£12 £15 64%^A (-27 + sVStt + 2t:'^)/SItt 

(3.33) 



with the remaining coefficients being determined by symmetry. The small fluctuation matrix is therefore 
given by 



Sttv aX 



( H 1 1 1 1 2-/X \ 

1 ^ 1 2-n 1 1 

1 1 /z 1 2-11 1 

1 2-n 1 /i 1 1 

1 1 2-/i 1 /i 1 

/ 



(3.34) 



\ 2-^JL 1111^ 

where ^ = 2(57r^ — 6-\/37r + 54)/97r^ ~ 1.59f9. (the A/a dependence has disappeared since the couplings a 
and A only appear in the dimensionless combination v^oA. The eigenvalues of the above matrix (neglecting 
the constants in front of it) are (with multiplicity 2), 2(/i — 1) (with multiplicity 3) and 6 (with multiplicity 
1). The zero modes correspond to flat directions for which deformations of the edge lengths leave the lattice 
geometry unchanged. The multiplicities of the eigenvalues are agree with the dimensions of the irreducible 
representations of the symmetry group of the tetrahedron. 




Fig. 10. Octahedral tessellation of the two-sphere, with arbitrary edge length assignments. 

For the octahedron (see Figure 10) one obtains instead the following coefficients of the small fluctuation 
matrix 

e?2 ^ 2VaA(216-12V37r + 57r^)/277r 
ei2 ei3 8\/aA (-27- 3\/37r + 27r^)/277r 

4\/aA (54 + 7r^)/97r 
£12 634 SVoA (-54 + 3V37r + 7r2)/277r 

ei2 £46 4\/oA (108 + 12VI-K + 7r2)/277r 

(3.35) 



again with the remaining coefficients being determined by symmetry. Up to a common factor of 2\fa\/27TT, 
the eigenvalues of the 12 x 12 small fluctuation matrix are given by 367r^ (with multiplicity 1), 972 (with 
multi " " ' ^ 




Fig. 11. Icosahedral tessellation of the two-sphere, with arbitrary edge length assignments. 

Finally, for the icosahedron (shown in Figure 11) one computes the following coefBcients of the small 
fluctuation matrix 



42 - 


16a 


/oA (270 - eVSTT + 7r2)/1357r 




16\ 


/aX (-675 - SOa/Stt + 87r^)/6757r 


ei2 ei4 - 


16\ 


/oA (270 - 6 a/3 + 7r^)/1357r 


£12 €34 - 


32a 


/oA (-675 + 15A/37r + 27r2)/6757r 


ei2 £45 - 


16a 


/oA (-675 + 15V37r + 27r2)/6757r 


£12 £38 - 


16a 


/oA (-675 + 15A/37r + 27r^)/6757r 


£12 £48 - 


16a 


/oA (675 + SOVStt + 7r^)/6757r 



(3.36) 

with the remaining coefficients being determined by symmetry. Up to a common factor of Sy/aX/GTb-K, the 
eigenvalues of the 30 x 30 small edge length fluctuation matrix are given (numerically) by 12340.173 (with 
multiplicity 3), 7238.984 (with multiplicity 5), 888.264 = 907r2 (with multiphcity 1), 20.887 (with multiplicity 
3), and zero (with multiplicity 18). 



The presence of the zero modes is related to the gauge (difFeomorphism) invariance of the gravitational 
action. The previous results can in fact be summarized as 

Tetrahedron (A^o 4) : 2 zero modes 
Octahedron(A*'o = 6) : 6 zero modes 
Icosahedron(A'o — 12) : 18 zero modes 

(3.37) 

If the number of zero modes for each triangulation of the sphere is denoted by Nz.m., then the results can 
be re-expressed as 

N^.rn. = 2No - 6 , (3.38) 

which agrees with the expectation that in the continuum limit, A^o ^ oo, N^.m./^Q should approach the 
constant value d in d space-time dimensions, which is the number of local parameters for a diffeomorphism. 
On the lattice the diffeomorphisms correspond to local deformations of the edge lengths about a vertex, 
which leave the local geometry physically unchanged, the latter being described by the values of local lattice 
operators corresponding to local volumes, and curvatures. The lesson is that the correct count of zero modes 
will in general only be recovered asymptotically for large triangulations, where A^o is roughly much larger 
than the number of neighbors to a point in d dimensions. It should be possible to find a similar pattern in 
higher dimensions. 

3.2 Edge Lengths as Metric Components 



Returning to the weak field expansion about flat space, it is easy to see that the above lattice gauge 
transformation corresponds to the diffeomorphisms in the continuum. Using the relationship between the 
metric perturbations and the edge length variations, obtained by choosing coordinates axes along the edges 
(as in Eq. ( |3.16 ) and Figure 7), 



6ll{n) 



5ll{n)) 



\{5ll{n)-5ll{n)-5ll{n)) 



sqin) 



one obtains from Eq. ( p. 23 ) the result 



Sgii = 5ll = 2xi(?^) - 2xi(ri + Ai) ~ -2^1X1 

^322 = Slj = 2x2(»^) - 2x2(»^ + A2) ~ -2 92X2 

Sgi2 = \{5ll- 511-611) 

= Xi{n + Ai) - Xi{n + Ai + A2) + X2{n + A2) - X2{'n + Ai + A2) 



(3.39) 



~ -^2X1 - dlX2 , 

(3.40) 

which can then be combined into the single famihar expression 

Sgtiu = -df.Xu - d^Xti , (3-41) 

and which is indeed the correct gauge variation in the weak field limit. 

Conversely, the above form of the lattice gauge transformation, Eq. ( 3.23| ), can be obtained from the 



form of infinitesimal difFeomorphisms in the continuum. In order to see this result, start from the definitions 
of diffeomorphisms 

9U^') = £r||^ffp-(=^) ' (3-42) 

as transformations which leave the infinitesimal line element, as well as any other coordinate invariant 
quantity, unchanged 

= g'^^{x')dx'^'dx"' = ds^ = gp„{x)dxPdx'' (3.43) 
under an arbitrary change of coordinates 

x'^ = x^' + x''{x) . (3.44) 

For infinitesimal variations one obtains 

g'tLui^') = 9iiu{x) + 6gf,^{x) = gf,,,{x) - Of^xix) d^x^{x) - gxu{x) df,x^{x) + O(x^) ■ (3.45) 

The above relationships express the well-known fact that metrics related by a coordinate transformation 
describe the same physical manifold. In the discrete case it reflects the invariancc of the lattice action under 
local deformations of the simplicial manifold which leave the local curvatures unchanged Since the 

continuum metric degrees of freedom correspond on the lattice to the values of edge lengths squared, one 
would expect to find analogous deformations of the edge lengths that leave the lattice geometry invariant, 
the latter being specified by the local lattice areas and curvatures, in accordance with the principle of 
discussing the geometric properties of the lattice theory in terms of lattice quantities only. Clearly the 
distance between lattice vertices will change under such a transformation, in accordance with the fact that 
only distances between fixed points will remain the same. This invariance is spoiled by the presence of the 
triangle inequalities, which places a constraint on how far the individual edge lengths can be deformed. In 
the perturbative, weak field expansion about a fixed background the triangle inequalities are not seen to any 
order in perturbation theory, they represent non-perturbative constraints. 



These considerations are further illustrated by the following elementary example pq l . In one dimension 
(zero space, one time dimension) one can discretize the line by introducing N points, joined by segments of 
lengths In- The only invariant term in one dimension is obviously the length of the curve, 

N 
n=l 

From the expression for the invariant line element, ds^ = gdx^ , one naturally associates g{x) with Z^, and 
the coordinate increment with the lattice spacing, dx ~ 1. One can take the view that distances can only 
be assigned between vertices which appear on some lattice in the ensemble, although this is not strictly 
necessary as distances can also be defined for locations that do not coincide with any specific vertex. 



The above one-dimensional action has an exact local invariance (compare with Eq. (3.23) in two dimen- 
sions) 

5ln = Xn+l - Xn , (3.47) 

where the Xn's represent continuous gauge transformations defined on the lattice vertices (actually, in order 
for the edge lengths to remain positive, one should also require Xn — Xn+i < In, which is satisfied for 
sufficiently small x's)- These transformations are in fact remarkably close in structure to the ones found 



in two dimensions (see Eq. (3.24)). Physically, the local invariance reflects the re-parametcrization, or 
coordinate invariance, of the original continuum action L ^ J dx\J g{x). It is the discrete form of the 
change Sg{x) — 2gdx- Variations of the edges which satisfy Eq. ( |3.47|) leave the physical length of the curve 
unchanged. In addition, given any two points on the curve, independent local gauge transformations can 
be performed on any of the vertices situated between the two points, while at the same time maintaining 
the same physical distance between them (which, for any assignments of edge lengths, is simply obtained by 
adding up the intervening edge lengths). It justifies the name lattice diffeomorphisms for the transformations 
of Eq. ( |3l^ ). 



In two dimensions one starts from the relationship of Eq. (3.16) between the squared edge lengths and 
the metric, and uses the expression for metric perturbations given in Eq. (3.39). From Eq. (3.45) one obtains 

5ll = Sgn = -2qdix' + {ll + ll-lDdix' 
Sll = Sg22 = -2qd2x' + ill + ll-ll)d2x' 
^^3 = '^.911 + ^.922 + 2(5512 

= i-ll + ll - ll) idix' + d2x') + {ll - ll - ll) {dix' + d2x') 

(3.48) 

^ There is here an (incomplete) analogy with ordin ary l attice gauge theories, in the sense that if one defines Un = e'" and 



Vn = e^"- , then the gauge transformation law of Eq. (3.47) can be re-written as Un — ► V^^UnVn+i, which parallels the gauge 
transformation law for the SU(N) gauge fields (/„^ in d dimensions, Un/j, — > V^^Unfi,Vn+fi., where are arbitrary N X N 
SU{N) matrices. The "update" Un ^ V ■ Un corresponds here to In ^ In + Sin- 



After introducing appropriate finite differences for the fields one then has 

511 = -2ll{xl-xl) + {ll + li-ll)ixl-xf>) 
5ll = -2ll{xl-xl) + {ll + ll-ll){xl-xl) 
5il = i^lj + il ll) ixl - xl) + HI ll ll) (X^ - X^) 

(3.49) 

Here upper indices on x label the components, while lower indices indicate the position. Taking Xi = X2 = 
X3 = 0; Xo 7^ 0, as well as h = h = 1, ^3 = (as appropriate for the square lattice), one finally obtains 
the simple result 

511 = 2x1 
5ll = 2x1 
511 = 2{xl + xl) 

(3.50) 



which is indeed equivalent to Eq. ( |3.23 ). This shows that the zero modes described in Eq. (3.23) correspond to 



lattice diffeomorphisms. (The term lattice coordinate transformations would appear to be equally suitable, 
provided one identifies the directions associated with the edges with a preferred coordinate system, and 
identifies changes in these coordinates as corresponding to variations in the squared edge lengths which leave 
the local curvatures unchanged). The case of flat space is obviously the simplest. By moving the location 
of the vertices around in flat space, one can find a different assignment of edge lengths which represents the 
same flat geometry. This leads to a d • iVo-parameter family of transformations for the edge lengths in flat 
space, and to a set of equivalent metrics which are all related by lattice diffeomorphisms, i.e. deformations of 
the edge lengths which leave the local curvature invariants unchanged, and respect the triangle inequalities. 

In conclusion, the previous analysis shows a direct correspondence between gauge transformations on the 
simplicial lattice 

q q + SlUx) , (3.51) 
and the analogous diffeomorphisms in the continuum 

g,,,y{x) — > 9ii.v{x) + 5gf,i,{x{x)) . (3.52) 

These transformations, in which suitable deformations of the edge lengths are shown to correspond to the 
local gauge transformations, should be contrasted with the set of what can be called trivial coordinate 
transformations. For a given assignment of edge lengths, introduce an arbitrary coordinate system, and a 



corresponding flat metric, within each triangle (or simplex in higher dimensions). Coordinate changes can 
then be performed within any triangle such that 

(3.53) 

within each triangle. These difTeomorphisms are trivial, in the sense that they correspond to a change in 
an arbitrary coordinate system, which was not part of the theory to begin with, as Regge's lattice theory 
is formulated exclusively in terms of coordinate- independent squared edge lengths, and not piecewise flat 
continuum metrics, which are highly degenerate. 

The confusion between the two types of invariance is, in our opinion, at the root of the erroneous 
conclusions drawn in Ref. , where it is argued that Regge gravity always needs a non-local gauge flxing 
term, to compensate for the fact that in the functional integral for gravity the integration is over "invariants" , 



the edge lengths. As the edge lengths correspond to components of the metric (see Eqs. (2.7), (3.18), (3.21) 



( |3.23 )), this cannot be true. The above discussion shows in detail that the situation is more subtle, and that 



a non-local additional, and in our opinion ad-hoc, gauge fixing term will most likely lead to an incorrect 
weighting, as already pointed out in ]40| . 

4 Arbitrary Curved Backgrounds 



The previous discussion dealt with the case of an expansion of the gravitational action about flat space, 
or a regular tessellation of the sphere, a manifold of constant curvature. To complete our discussion, we now 
turn to the more complex task of exhibiting explicitly the local gauge invariance of the simplicial theory, for 
an arbitrary background simplicial complex. To this end we write 

I' = ll^+q^+Sl', , (4.1) 

where qi describes an arbitrary but small deviation from a regular lattice, and Slf is a gauge fluctuation, 
whose form needs to be determined. We shall keep terms 0{q^) and 0{q SP), but neglect terms 0{5l*). 

The squared volumes V^^ (cr) of n-dimensional simplices a are given by homogeneous polynomials of order 
(Z^)". In particular for the area of a triangle with arbitrary edges I1J2, h one has 

SAl = li-q + q + ll) 5ll + - ll + ll) 5ll + + ll - ll) 5ll , (4.2) 



and similarly for the other quantities in Eqs. (2.16), (2.17) and (2.18) which are needed in order to construct 
the action. 

For our notation in two dimensions we will refer to Figure 12. The subsequent Figures 13 and 14 illustrate 
the difference between a gauge deformation of the surface which leaves the area and curvature at the point 
labeled by invariant, and a physical deformation which corresponds to a re-assignment of edge lengths 
meeting at the vertex such that it alters the area and curvature at 0. In the following we will characterize 
unambiguously what we mean by the two different operations. 




Fig. 12. Notation for an arbitrary simplicial lattice, where the edge lengths meeting at the vertex have been deformed 
away from a regular lattice by a small amount qi (m,inim,ally deformed equilateral lattice). 




Fig. 13. Local gauge deformations of the lattice act on the edge lengths meeting at the vertex 0, and are performed in 
such a way that the area and curvature at the vertex are left unchanged. 



Fig. 14- Physical deformations change the area and curvature at the vertex 0, thus changing the lattice geometry. 



4.1 Equilateral Lattices 

First consider the expansion about a deformed equilateral lattice, for which loi = 1 to start with. A 
motivation for this choice is provided by the fact that in the numerical studies of two-dimensional gravity the 
averages of the squared edge lengths in the three principal directions turn out to be equal, {if) = {I2) = {l^)- 
The baricentric area associated with vertex is then given by 

4^01 + 902 + 916 + 912) 

4902 + 903 + 912 + 923) 

4903 + 904 + 923 + 934) 

4904 + 905 + 934 + 945) 

4905 + 906 + 945 + 956) 

4906 + 901 + 956 + 916) 

(4.3) 

Our normalization is such that Aq = ^ for = 0. Equivalently one can write, in more compact notation, 
at the vertex 

A = Ao{q) + ^va{<i) ■ + 0{,5f ) , (4.4) 
with 8P = {SIqi, . . . , 61qq). After adding the contributions from the neighboring vertices one obtains 

A= J2 Mq)+ Mq) ■ + 0{5f) . (4.5) 
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0{6l^) 



Therefore the area associated with the vertex wih remain unchanged provided the variations in the squared 
edge lengths meeting at satisfy the constraint 



VAiq) -SP ^ 



(4.6) 



This is nothing but the curved space equivalent of the condition of Eq. (3.30), which for the flat equilateral 
lattice takes the form 

1=1 

Alternatively, if one considers a dual subdivision, one has to consider the dual area associated with vertex 
0. In this case one has 



A = Aoiq) 
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4 • 33/2 



(5/^1 (1 + 2qoe - 4goi + 2qo2) 



(1 + 2goi - 4(702 + 2(703) 
(1 + 2go2 - 4go3 + 2^04) 



+ Sll^ (1 + 2(/o3 - 4(/o4 + 2(705) 
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(1 + 2go4 - 4go5 + '^qoe) 
(1 + 2go5 ~ 4go6 + 2(701 ) 



(4.8) 



leading to a result formally similar (in fact in this case identical) to the baricentric case. 

A similar calculation can be done for the curvature associated with vertex 0. One has for the deficit 
angle at 

(5 = (5o(g) + ^-^ [ (5/gi (3-2(7o6-goi-2go2 + gi6 + 9i2) 



+ (5^02 (3 - 2(7oi - qo2 - 2go3 + '7i2 + to) 



+ 



(4.9) 



and therefore for the variation of the sum of the deficit angles surrounding 



and 



^(e^'^I = e 

\ h I Po - Pa 



E E So{q) + vn{q)-SP + 0{6l') 

Po...Pe Po-Pe 



(4.10) 



(4.11) 



with in this expected, VR{q) = 0. 

FinaUy for the curvature squared associated with vertex one computes 

4 r 

6'^ /A = So/Aoiq) + Sl^, (goi + 902 + qo3 + qo4 + 905 + 906 - 9l2 - 923 - 934 - 945 - 956 " 916) 

+ ^^02 (901 + 902 + 903 + 904 + 905 + 906 - 912 - 923 - 934 - 945 - 956 - 9l6) 



0{SI^) 



Adding up aU seven contributions one gets 



and therefore 



\ h ) Po-Pe 



E SyA= E {6'/A), + vn.iq)-Sp + Oi5l'') 

Po-.-Pe Po-.Pe 



(4.12) 



(4.13) 



(4.14) 



In this case the curvature squared associated with the vertex will remain unchanged, provided the variations 
in the squared edge lengths meeting at satisfy the constraint 



(4.15) 



which provides a second constraint on the edge length variations SP at the vertex 0. Again this constraint 
is the generalization to curved space of the condition of Eq. ( 3.26| ) , which was valid for flat space. 



4.2 Square Lattice 



A similar calculation can be performed for the square lattice with ^oi = ^02 = 1 and /03 — V^- The 
baricentric area associated with vertex is now given by 
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(4.16) 



P0...P6 Po-.-Pe 



Our normalization here is such that Aq = 1 for qi = 0. Summing up all relevant contributions, one can write 

(4.17) 

leading to the invariance constraint on the edge length variations 

v'A{q)-SP = . (4.18) 

This provides a first constraint on the edge length variations 6P at the vertex 0, for the deformed square 
lattice. 

For the dual (Voronoi) area associated with vertex one has instead 
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(4.19) 



leading to a constraint similar to the one for the baricentric area. For the curvature associated with vertex 

one has 

S = 5o{q) + \ 
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(4.20) 



which gives now 



E ^= E So{q)+v'n{q)-SP + 0{5l^) , 

Po.-.Pe Po-Pe 



(4.21) 



with v' R{q) = 0. 

Finally, the curvature squared associated with just the vertex is given by 

g^/A = (5o/^o(g) + \ 5^01 (901 + 902 + 903 + 904 + 905 + 906 - 912 - 923 - 934 - 945 - 956 - 916) 
+ <^^02 (901 + 902 + «03 + 904 + 905 + 906 - 912 - «23 " 934 - «45 " 956 " 9l6) 



+ 



+ 0{5l^) 



(4.22) 



The labeling in the previous formulae is a bit clumsy in dealing with nearest next-nearest neighbor interac- 
tions; for the labeling in the following formula we refer to Figure 15. 
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Fig. 15. Labeling of lattice vertices for the expansion of the curvature squared term around a minimally deformed 
square lattice. 

Adding up the contributions from the seven distinct vertices — 6 one obtains for the term linear in 51"^ 



<^^o,i( 4go,i + 2g'o,2 + 90,-1 - 9o,-2 - 9o,3 + 2^0,-3 

-91,3 - 2g2,3 - 9-1,-3 - 2g_2,-3 - 9-i,2 - 9i,-2 
+91,1+1 + 2gi4_2 + 2qi_i+3 — (7i+i_i_2 — qi+i^i+z 

+92,2-1-3 — 2g_2,l-2 — 9-2,-2-2 — 9-2,-2-3 + 9-2-2,-2+1 + 9-2-2,-2-3 

— 2g3^1+3 — 93,2+3 — 93,3+3 — 9-3,-2-3 + 93+3,1+3 + 93+3,2+3) 

+ ^^o,3( -90,1 - 90,2 + 4^0,3 + 2gro-i + 2g'o,-2 + 90,-3 

-91,3 - 92,3 - 9-1,-3 - 9-2,-3 - 29-1,2 - 291, _2 

— 91,1+1 — 9l,l-2 — 291,1-1-3 + 9-1,2-1 + 91-1-1,1-2 + 91-1-1,1-1-3 

— 92,2+2 — 92,2-1 — 292,2+3 + 9-2,1-2 + 92+2,2-1 + 92+2,2+3 



+2g3,i+3 + 2g3,2+3 + 93,3+3 — 93+3,1+3 — 93+3,2+3) 



+ 



(4.23) 



which can be written in short form as 



(4.24) 
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Again the curvature squared associated with the vertex will remain unchanged provided the variations in 
the squared edge lengths meeting at satisfy the constraint 



which provides the second constraint on the edge length variations SP at the vertex 0, for the deformed 



In conclusion, we have shown explicitly how gauge variations of the edge lengths at each vertex can be 
defined by requiring that the action contributions be locally invariant. We have looked at small deformations, 
but large deformations can be treated as well along the same lines, provided one is careful not to violate the 
triangle inequalities, which impose a non-pcrturbativc cutoff in orbit space The same approach can also be 
extended to higher dimensions, leading to similar (but rather more complicated, when written out explicitly!) 
results. The main conclusions do not change. 



In the previous section we have discussed the invariance properties of the lattice action for pure gravity. 
Next a scalar field is introduced, as the simplest type of dynamical matter that can be coupled to gravity. 
The scalar lattice action in the continuum is 



The dimensionless coupling ^ is arbitrary; two special cases are the minimal (^ = 0) and the conformal 



w'fl2 (q) -SP = . 



(4.25) 




5 Scalar Field 




(5.1) 



{£_ — i) coupling case. In the following we will mostly consider the case ^ = 0. 



5.1 Construction of the Lattice Action 



On the lattice consider a scalar (j>i and define this field at the vertices of the simplices. The corresponding 
lattice action can be obtained through the usual procedure which replaces the original continuum metric 



with the induced metric on the lattice, written in terms of the edge lengths g|. Here we shall consider 
only the two-dimensional case; the generalization to higher dimensions is straightforward. It is convenient 
to use the notation of Figure 16, which will bring out more readily the symmetries of the resulting lattice 
action. Here coordinates will be picked in each triangle along the (1,2) and (1,3) directions. 




Fig. 16. Labeling of edges and fields for the construction of the scalar field action. 

To construct the scalar lattice action, one performs in two dimensions the replacement 

,,.(.) ,.,(A) ^ ( , f ^ f + ) , (5.2) 

which then gives 

deig^,{x) det5,,(A) = \ {2{llll + llll + llll) -lt-ll-lt] = \A\ , (5.3) 

and also 

'T^^ - •<'*'<^) - W ( \ - '■!) "I " ® ) ^ 

For the scalar field derivatives one writes |4^, 

which corresponds to introducing finite lattice differences defined in the usual way by 

d^(f> — > (A^(/))j = (j>i+^ - (j>i . (5.6) 

Here the index fi labels the possible directions in which one can move from a point in a given triangle. The 
discrete scalar field action then takes the form 



1^(1^ ^) ^ TeT.J^ ['i(<^2 - ^iKh - 0l) + ll{^3 - ^2K4>1 - 02) + lli^l - 03)(<^2 - 03)] • (5.7) 

Using the identity 

(0. - 0,)(0, - 0fe) = i [{<P^ - <^jf + (0, - <Pkf - (0, - <t>kf] , (5.8) 



one obtains after some re-arrangements the simpler expression 

I^Al'A) = ^E^^^l^)' ' (5-9) 
<ij> 

where Aij is the dual (Voronoi) area associated with the edge ij. In terms of the edge length lij and the 
dual edge length hij, connecting neighboring vertices in the dual lattice, one has Aij = \hijlij (see Figure 
17). Other choices for the lattice subdivision will lead to a similar formula for the lattice action, but with 
the Voronoi dual volumes replaced by their appropriate counterparts in the new lattice. 




Fig, 17. Dual area associated with the edge li (shaded area), and the corresponding dual link hi. 
For the edge of length li the dihedral dual volume contribution is given by 



16^123 16^234 



with hi is the length of the edge dual to li. The baricentric dihedral volume for the same edge would be 
simply 

Al, = (Ai23 + v4234)/3 . (5.11) 

It is well known that one of the disadvantages of the Voronoi construction is the lack of positivity of the dual 
volumes, as already pointed out in 0|. Thus some of the weights appearing in Eq. (|]^) can be negative for 
such an action. On the other hand, for the baricentric subdivision this problem does not arise, as the areas 
Aij are always positive due to the enforcement of the triangle inequalities. It is immediate to generalize the 
action of Eq. (|5.9| ) to higher dimensions, with the two-dimensional Voronoi volumes replaced by their higher 
dimensional analogs. 

Mass and curvature terms can be added to the action, so that the total scalar action contribution becomes 

< ij > i 



The term containing the discrete analog of the scalar curvature involves the quantity 

AR, = Y,6h - ■ (5.13) 

In the above expression for the scalar action, Aij is the area associated with the edge lij , while Ai is associated 
with the site i. Again there is more than one way to define the volume element Ai, but under reasonable 
assumptions, such as positivity, one expects to get equivalent results in the lattice continuum limit. In the 
following we shall mainly consider the simplest form for the scalar action, with = ^ = 0. 

One of the simplest problems which can be studied analytically in the continuum as well as on the lattice 
is the analysis of small fluctuations about some classical background solution. In the continuum, the weak 
field expansion is often performed by expanding the metric and the action about flat Euclidean space 

5^"^ ^ S^^ + K h^^ . (5.14) 

In four dimensions n = V327rG, which shows that the weak field expansion there corresponds to an expansion 
in powers of G. In two dimensions this is no longer the case and the relation between k and G is lost; instead 
one should regard k as a dimensionless expansion parameter which is eventually set to one, k = 1, at the end 
of the calculation. The procedure will be sensible as long as wildly fluctuating geometries are not important 
in two dimensions (on the lattice or in the continuum). The influence of the latter conflgurations can only 
be studied by numerical simulations of the full path integral |jri|, ^ . 

In the continumn, the Feynman rules are obtained by expanding out the action in the weak fields h^^{x), 

\ J d^x^gf^^d^q^d^cl) ^ \ j d^x{d<pf + \ J d^xh^,{^S^,{d<l>f~d^cl)d,^}+0{h^) , (5.15) 

and by then transforming the resulting expressions to momentum space. 

On the lattice the action is again expanded in the small fluctuation flelds , which depend on the specific 
choice of parameterization for the flat background lattice - a convenient starting point is (in two dimensions) 
the square lattice with diagonals. It is convenient to define the edge variables at the midpoints of the links 
p8[ . For the edge lengths one then defines the lattice Fourier transforms as 

= ££(^ 

- £ £(S^ ^='*>- 

(5.16) 

while the scalar fields are still defined on the vertices, and are Fourier transformed in the usual way, namely 



These formulae are completely analogous to the ones used in developing the perturbative expansion for 



lattice gauge theories |44 . They are easy to generalize to higher dimensions when a simplicial subdivision 
of a hypercubic lattice is employed, as first suggested in Following this procedure, the Feynman rules 
for the lattice scalar field action were derived in JTsf (using the Voronoi dual volumes) , and shown to agree 
completely with the continuum Feynman rules for small momenta. For more details, and the computation of 
the Feynman diagrams relevant to the conformal anomaly, we refer the interested reader to the cited work. 

5.2 Equilateral Lattice 



In order to compare and analyze the difference between the two volume discretizations (baricentric vs. 
dual) for the scalar action, it will be useful to look at their form for small deformations of the edges about 
a regular lattice, such as an equilateral or a square one. One would expect that the precise form of the 
discretization should not matter, as long as the correct long distance (small momentum) properties are 
preserved. Let us see how this can come about. 



The next step is therefore the expansion of the lattice scalar field action of Eq. (5J) with volumes Aij 
defined via for example a baricentric subdivision, starting from an equilateral lattice with If = l^'^ + Slf, 
and 1^ = 1. It will be sufficient to limit oneself to the contributions coming from one site (0) and its six 
neighbors (1-6), which for the kinetic term is given by 
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(5.18) 

For our notation and labeling of the edge lengths we refer again to Figure 12. In the case of the dual 
(Voronoi) subdivision, the results is identical to this order except for the replacement of the coefficient 
2471 ~* 12^/3' '^^^ interpreted as a rescaling of the gravitational coupling between the metric field 

and the scalar. For fields that are smoothly varying on the scale of the cutoff, the scalar action contribution 
is invariant under a gauge transformation on the edges acting at the origin 0, if the defining condition for 
gauge transformations, analogous to Eq. ( |3.26 ) pertaining to the square lattice in the weak field expansion, 



is satisfied 

6 

J2 Slfin) « . (5.19) 

This result is also very similar to what happens in one dimensions, where the exact invariance of the scalar 
action can be written down explicitly with more ease p8| . 

For the same (equilateral) lattice let us consider next the mass term; it is given by 



I h4 



(5.20) 



Again expanding about the equilateral lattice with If = /° ^ + dtf one obtains for the baricentric subdivision 
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(5.22) 



2 12n/3 

Again to this order the mass term is invariant under a gauge transformation on the edges acting at the 
origin 0, if the defining condition for gauge transformations (compare to Eq. ( |3.30 ) obtained in the weak 
field expansion) is satisfied 

6 

i=l 

Finally one can consider the curvature term, 



Its expansion is given by 
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(5.24) 



(5.25) 



with baricentric and dual forms identical to all orders in since Ai does not appear in this term. Again 
it is obvious that this term is invariant under gauge transformations at 0. 
Further higher order terms involving the curvature, such as 

A., 



(5.26) 



become increasingly complicated in their structure, as they involve neighbors which are further apart. They 
are strongly suppressed for smooth manifolds due to the presence of the deficit angle. In the baricentric case 
one finds 



M'5^( 2< + 0-5/o'2 



• 5ll^ ~ 6ll, - Sll^ - Sll, - • 6ll, + •••) + 0{SI') 



(5.27) 



and we have omitted the terms involving (02 — (f'o)'^ etc. since they can be obtained by symmetry. Again 
for slowly varying scalar fields the above term is invariant under gauge transformations at the origin. 

5.3 Square Lattice 



A similar calculation can be performed for the square background lattice with again If = Z?^ + 5lf, and 
^01 = ^02 = 1 and ^03 = \/2- Using the baricentric subdivision, the kinetic term gives 
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(5.28) 

Again for smoothly varying scalar fields, the above action will be invariant under gauge variations of the 
edge lengths at the vertex 0, provided one has 61qi + 2(5Zq2 + 61^^ + SIq^ + 261^^ + 61qq — (which is the scalar 
field action analog of Eq. ( 4.18| )). 

In the case of the dual subdivision, one obtains for the same lattice 
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(5.29) 

Due to the asymmetry of the lattice in this case, the difference between the two discretizations is quite 
marked. We will argue below that the dual volume form represents in fact an "improved" discretization over 
the baricentric form. 

For the mass term one proceeds in the same way, and for the baricentric subdivision one obtains 
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(5.30) 



Again, this term is invariant under a gauge transformation if Eq. ( 4.1S| ) is satisfied. For the dual (Voronoi) 
case one finds similarly 
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which is invariant again if the area at is invariant 
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Finally for the curvature term of Eq. ( 5.24 ) one has 
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(5.32) 



with the baricentric discretization equal to the dual discretization to all orders in SP, since the local area 
Ai does not appear. Again it is obvious that this term is invariant under gauge transformations at 0. 

Let us now return to the apparent discrepancy between the scalar kinetic term in the baricentric and 
dual discretizations. It is useful to look at the two discretized forms in momentum space. Thus we assume 
that the edge length fluctuations = 51^/21^"^ and cj) at the point i, j steps in one coordinate direction and 
k steps in the other coordinate direction from the origin, are related to the corresponding and (p at the 
origin by 
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where LOi = e , and ki is the momentum in the direction i. Similarly for </> one writes 
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where to'i — 



After making the substitution — s- of Eq. ( 3.1j ) (which is the same as the change of variables in 



Eq. (3.9)) one obtains for the interaction in momentum space, in the dual case 



pI'P'' [l + ji~hn + h22) + 0{h^) 
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(5.35) 



which indeed coincides with the weak field expansion of the kinetic term for the original scalar field action 



in the continuum, Eqs. (xl) and ( 5.15 ). In the baricentric case one obtains instead 

pI^^ [h + + ^^22 - hl2) + 0{h^) ] 



+ pI^"" [5 + g( + /lll-^22-/il2)+0(/l2)] 



(5.36) 

which at first looks quite different from the dual (Voronoi) volume case. 




Fig. 18. Two equivalent triangulation of flat space, based on different subdivisions of the square lattice. 

On the other hand this result is not too surprising, as an analogous situation occurs in flat space. Consider 

for example the two lattices of Figure 18 (up to now we only considered the one corresponding to 18a), and 
compute in each case the action for one momentum mode, which is just the inverse scalar propagator in 
momentum space. In the absence of any interaction terms along the lattice diagonals one has 

2 ^ [ 1 - cosp^ ] ^ pI+pI + 0(/) . (5.37) 
ft 

The interaction terms along the lattice diagonals in Figure 18a contribute 

2[ 1 - cos(pi +P2)] pI+pI + 2piP2 + 0(p^) , (5.38) 
while the interaction terms along the lattice diagonals in Figure 18b contribute 

2[ 1 - cos(}Ji -P2)] ~ pI+pI- 2pip2 + Oip^) . (5.39) 

Thus when one averages over the two equivalent contributions (since the lattice is dynamical, and both 
contributions are equally probable) one obtains 

2[l-cospicosp2] ~P?+P2+0(/) , (5.40) 



which is now rotationaUy invariant to order p^. When the same procedure is apphed to the lattice scalar 
action in the presence of the gravitational field, Eqs. ( 5.35D and ( 5.36| ) can be shown to become equivalent, 



after a rescaling of the gravitational coupling. Still, the action based on dual Voronoi volumes of Eq. ( 5.35 ) 
appears to lead to a more attractive discretization, as the unwanted terms do not appear at all for the choice 
of lattice of Figure 18a, and no averaging over the gravitational field fluctuations is necessary to exhibit the 
correct correspondence with the continuum action. How seriously one takes this class of problems depends 
on how seriously one trusts low-order perturbation theory about a fixed flat background as a tool for the 
study of fluctuating geometries. 

5.4 Invariance Properties of the Scalar Action 



It is instructive to look at the invariance properties of the scalar action under the continuous lattice gauge 



transformations defined in Eqs. ( p.23| ) and ( 3.49 ). Physically, these local gauge transformations, which act on 
the vertices, correspond to re-assignments of edge lengths which leave the distance between two fixed points 
unchanged. In the simplest case, only two neighboring edge lengths are changed, leaving the total distance 
between the end points unchanged. On physical grounds one would like to maintain such an invariance also 
in the case of coupling to matter, just as is done in the continuum. 

The scalar nature of the field requires that in the continuum under a change of coordinates a; — > x', 

</)'(x') = 0(x) , (5.41) 

where x and x' refer to the same physical point in the two coordinate systems. Let us first look at the 
one-dimensional case, which is the simplest. On the lattice, as discussed previously, gauge transformations 
move the points around, and at the same vertex labeled by n we expect 

0„ -> 0; « 0„ + f 'tt±lfl±^\ e„ , (5.42) 



One can determine the exact form of the change needed in 0„ by requiring that the local variation of the 
scalar field action 

— (0,1 + A0„ - 0„-l)^ + — — {(pn+l ~ 4>n- ^<t>nf 

-j^ (0„ - <^«-i)' - ^ (0n+i - </)«)' = (5.43) 

be identically zero. Solving the resulting quadratic equation for A(/)„ one obtains a rather unwieldy expres- 
sion, which to lowest order is given by | ]38| 



and which is indeed of the expected form (as well as being symmetric in the vertices n — 1 and n + 1). For 
fields which are reasonably smooth, this correction is suppressed if It/iri+i — (^nl/ln <C 1. On the other hand it 
should be clear that the functional measure dcfin is no longer manifestly invariant, due to the rather involved 
transformation property of the scalar field. 

A similar line of argument can be pursued in higher dimensions. Thus in two dimensions one should 
require that locally the variation of the action contribution be again zero for edge length deformations Shj , 
associated with edges meeting at the vertex i, and which correspond to lattice gauge transformations (in the 



weak field, for example, they have the explicit form given in Eqs. (3.21) to (3.23)), namely 




2 



where j labels the neighbors of the site i. The transformation law for (pi is then determined by solving the 
above equation for A(j)i, given an arbitrary gauge variation SP{x) at the vertex i. 

5.5 Equations of Motion and Lattice Energy Momentum Tensor 



The equations of motion for are obtained from 

d 



\^<ij> i 



= 



and read 



(/J — 



lij lij 



+ m'^Ai (pi = 



(5.46) 



(5.47) 



where the notation indicates that the site j is taken to be adjacent to i. For the choice of indices in 
Figure 12 (see also Figure 17), the equation with i = reads 
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(5.48) 
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and represents the discrete analog of df^g'^'^'^ gdiy(f> + m^(f> = 0. Eq. (5.47) can be re-written, for = 0, as 
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(5.49) 



The d can be interpreted as normalized hopping amplitudes, with a normalization 

condition 



It is easy to verify that this last property does not depend on the way the original lattice has been subdivided 



in order to construct the dual lattice. Both the Voronoi construction of Eq. (2.17) and the bariccntric one 



of Eq. ( 2.18| ) lead to the same normalization for . And indeed it is easy to show that this result holds in 



higher dimensions as well. 

Some very interesting properties regarding the spectrum of the Laplacian on flat random lattices have 
been worked out in |l^ and the reader is referred to these papers for further details. Perhaps one of the 
most interesting results which can be derived by the method of replicas in the weak disorder limit is that, at 
least in sufficiently low dimension, the spectrum of the Laplacian coincides with the continuum result at low 
frequencies, with subleading corrections which then reflect specific aspects of the edge length distribution 
such as its Poissonian form |24[ In other words, quenched random lattices and regular lattices give the 
same (continuum) result at low frequencies. 

In the continuum the energy momentum tensor for matter described by action is defined via the 
relationship 

SI,n = \ j d'^x^T'^'^Sg^, (5.51) 
For infinitesimal gauge variations, which have the form 

59t,y{x) = -g^,\{x)duX^{x) - gxu{x)df,x^{x) - dxgi_i,,{x)x^{x) = g^i. , (5.52) 



one expects 5I„i = 0. After an integration by parts in Eq. (5.51), one then obtains 



(n);^ = . (5.53) 



The energy momentum tensor defined by equation Eq. ( ^.51 ) will be conserved if and only if the matter 



action is a scalar. More concretely, for one real scalar field of mass m, one has 

T^. = c)M'/'9.0-i(7^.(5A<^aV + m202) (5 54) 

Taking the trace one obtains 

n = (^) ^^'/'S^'/'-i^v . (5.55) 

It is natural to proceed on the lattice in analogy with the continuum case. For the scalar action contribution 
of Eq. ( |5.9| ) one computes the variation 

= E E - ^^)' W (^] Sll ■ (5.56) 

k ij ^ \ / 




Fig. 19. Notation for Eq. ( 5.56 ) describing the variation of the action. 

Since the derivative term inside the sum is non-zero only for edges ij which belong to triangles touching 
the edge k (see Figure 19), only four terms contribute to the sum over < ij >. It is actually more convenient 
to start from the equivalent form for the scalar action given in Eq. (p^, and one obtains 
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(5.57) 
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Using the definition for the dual volumes (see Eq. ( 2.17 ), and Figure 20 for our notation here) one can 
re-write the above expression more compactly as 
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Therefore one can introduce the quantities such that 



(5.58) 
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(5.60) 



associated with the edge labeled by 1 in Figure 19, and similarly for all the other edges in the lattice. It is 
clear that this equation defines the analog of the energy- momentum tensor in the discrete case. 



'3 



Fig. 20. Labeling of the dual areas appearing in Eq. (5.5t ). 

The coefficients in tlie above expansion in Sl^ can be regarded as tlie components of the Regge lattice 
analog of the energy-momentum tensor, 

Tfe = T^. II , (5.61) 
just as we can define for the simplicial components of the metric tensor p5| 

gk = g^Xj-^ = If . (5.62) 



From Eq. (2.7) and Eq. (3.3£) 



sii + sil^-sq^ 



(5.63) 



and therefore within triangle 1, (with vertices 1,2,3 and choosing coordinates along 23 and 21), one has (see 
Figure 20) 

= Sgii 

= (5.922 

= Sgii + 5g22 - 2(5gi2 . 
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For one triangle (for example, triangle 1 in Figure 19) one obtains 



(5.64) 



(5.65) 



In the lattice case it is clear that, inserting for the variation of the squared edge lengths corresponding to 



gauge variations, as in Eq. (3.21) and Eq. (3.23), and then equating the resulting coefficients of the arbitrary 
gauge parameters xli to zero, gives the discrete equation of conservation for the energy-momentum tensor. 
On the other hand, the equations of motion for If (as opposed to the equations of motion for 0, which 



are given in Eq. ( 5.48 )) are obtained directly from Eq. ( 5.58 ), namely 
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It corresponds to the continuum equation = 0. In the presence of a cosmological constant term, A At^ 
an additional term appears in the equations of motion, which become, again with the notation of Figures 19 
and 20, 
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and relates the squared edge lengths to the derivatives (or, more properly here, finite differences) of the 
scalar field. 

6 Gravitational Functional Measure 



In this section we will re-examine the issue of the gravitational functional measure in light of the results 
of the previous sections, and in particular the local gauge invariance of the lattice gravitational action. It 
is well know that in ordinary (lattice) gauge theories the invariance of the action selects a unique measure 
(the group invariant Haar measure) . A natural way to construct the gravitational functional measure in the 
continuum is to introduce a metric over metrics (or super- metric), and then compute the resulting functional 
volume element. We shall see that what appears at first as a rather straightforward procedure, is in fact 
affected by a number of rather subtle ambiguities. 

6.1 Continuum Case 

Following DeWitt one introduces a super- metric G over metric deformations 5g^u{x)^ which in the 

simplest local form leads to the following norm-squared for metric deformations 

\\5gf ^ j d^xG^''^'^%{x)]5g^,{x)5g^p{x) , (6.1) 

with the inverse of the DeWitt supermetric given by 

G^'''^%{x)] = \^)[g^^%x)g'^\x)+g^^P{x)g''^{x) + \g^^^{x)g"P{x)\ , (6.2) 

and A ^ — 2/(i, to avoid the vanishing of the determinant of G. It is easy to check that the above expression 



for 1 1 (5.9 IP is invariant under difFeomorphisms. ^ The usual procedure is then to derive the functional measure 
in the form 

/ diJi[g] = j n(det[G(.9(x))])5 , (6.3) 

X fj,>U 

with the determinant of the super-metric G'-'''^''"'^ {g{x)) given by 

det G{g{x)) cx (1 + ^dX) [g{x)]'^'^-^^'^''+^^^^ . (6.4) 

Up to irrelevant constants, in four dimensions it reduces to the very simple expression 

/ d^l[g] = JH [g(x)]('^-4)('^+i)/« n d9,A^) ,"^4 / 11 11 '^9,A^) " (6-5) 

Unfortunately this measure is not gauge invariant, if the product over x is interpreted as one over 'physical' 
points, and coordinate invariance is imposed at one and the same 'physical' point, as discussed in p8[ . Here 
this is seen as a consequence of the fact that ||(5.g|p has been split in two separately non-invariant parts. The 
measure that does satisfies the invariance property is 

j = / n (det[G(5(x))/y^])' ' (6-6) 

and was originally proposed by Misner |Q. Explicitly, 

/ d^i[g] = I Yl n '^a.A^^ ■ (6.7) 

X 

It has the property of being scale invariant in any dimension. Unfortunately it is also singular, and needs 
therefore to be regulated in some way for small g (on the lattice this requires some cutoff for small local 
volumes) . 

Indeed both measures can be obtained as particular cases if one introduces a real parameter oj, and writes 
11-5511' = J d''x{g{x)r'^ G^''-"^[g(x);a;]%.(x)<55a/3(x) , (6.8) 

with 

G^^^"'^b(a;);c.] - \ {g{x)f-''^'^ [g^^%x)g'^0{x) + g^l'{x)g'^'^{x) + \g^^'^{x)g'^P{x)\ . (6.9) 

The metric in function space is obviously left unchanged by this rewriting, but the measure (obtained from 
det G) depends on u, and by the usual argument one obtains the parameterized functional gravitational 



measure 



/ ^^^[5] = /n [vsM]" n ■ (6.10) 

X ii>i/ 



^While it should be clear that the functional norm should be invariant, it is less obvious that it should be local. Ultimately 
the justification for locality lies in the fact that the resulting functional measure is local in the continuum, a rather desirable 
feature. 



with the measure parameter a related to the choice of lo via 

o = + . (6.11) 

For uj = one obtains the DeWitt measure of Eq. (^^), while for a; = 1 one has the Misner measure of 
Eq. ( |6.7[ ) . The close relationship between the DeWitt and Misner measures was pointed out in |Q . ^ 

As there is no way of deciding between these two choices, or any intermediate one for that matter, one 
is forced to consider a as an arbitrary (and hopefully ultimately irrelevant) real parameter of the theory, the 
only constraint being a > —{d + 1). In general the volume factor in the functional measure is absent 
in d dimensions for the special choice lu = 1 + A/d. It should be emphasized here that gauge invariance 
does not select a specific value for a, but does otherwise completely constrain the form of the measure. The 
criteria of simplicity and universality would suggest the above to be the preferred choice for lo. The reason 
for the ambiguity in the gravitational functional measure appears to be a lack of a clear definition of what 
is meant by J^^ . In spite of some recurrent claims to the contrary, such an ambiguity persists in all lattice 
formulations. Also, the volume term in the measure is completely local since it contains no derivatives. It 
does not effect the propagation of gravitons, as it contributes 6'^{Q) terms to the effective action. To some 
extent these can be regarded as renormalizations of the cosmological constant, since they clearly only affect 
the distribution of local volumes. As such they are expected to only affect the short distance behavior of the 
theory, leaving the more interesting universal large distance properties unmodified (the recent work in Ref. 
I^lf gives a more concrete realization of these ideas in the framework of continuum perturbation theory). 



Support to this interpretation also comes from a number of simple examples |38 
6.2 Lattice Transcription 



Let us examine the consequences of this discussion in the discrete case. As the edge lengths play the role 
of the metric in the continuum, one expects the discrete measure to involve an integration over edge lengths 
1^, 1^] . Indeed the induced metric at a simplex is related to the edge lengths squared within that simplex, via 
the expression for the invariant line element ds^ — g^j^^dx^dx^ . The relation between metric perturbations 
and squared edge length variations for a given simplex in d dimensions is (see Eqs. (2.7) and ( 3.39| ), and 
Figure 3) 

59r^) = \ {Sll, + < - Sl^) . (6.13) 

''The gravitational measure has to be modified in the presence of the matter fields. For an n-component massless scalar field 
the invariant measure is 

f rf/^M = f Ylg'^^H^) f[dr{x) , (6.12) 

X a — l 

and therefore a = -(d + 1) + - 1) ii^+ll + n/2. 



Thus for one simplex the integration over the metric is equivalent to an integration over the edge lengths, 
and one has 

, , CT d(d-l) d{d+l)/2 

(^vW^j U'^dvi^) ^ i-k)"^ [V'iiix n ■ (6.14) 

^ ' ^ i>j fc=l 

There are d(d + l)/2 edges for each simplex, just as there are d{d + l)/2 independent components for the 
metric tensor in d dimensions. (We are ignoring for the moment the triangle inequality constraints, which 
further require all sub-determinants of gij to be positive as well, including the obvious restriction If. > 0). 
The extension to many simplices glued together at their common faces is then immediate, and after summing 
over all simplices one obtains, up to an irrelevant numerical constant, 

/ df^r] = r II [V4s)r n < ■ (6-15) 

In four dimensions the lattice DeWitt measure (cr = 0) is particularly simple, 

J d^if] = ^"n<-^[^?>] • (6.16) 



Here F^[l] is a (step) function of the edge lengths, with the property that it is equal to one whenever the 
triangle inequalities and their higher dimensional analogs are satisfied, and zero otherwise. ^ 

The above lattice measure over edge lengths has recently been used extensively in numerical simulations 



of simplicial gravity |15| |5^, ^ |10|, ^ (Figure 21 gives an example of the edge length distribution 
obtained in four dimensions with measure parameter a — 0). One would expect that physical properties of 
the theory should not depend on the arbitrary parameter cr, but for the moment this universality argument 
remains largely an unproven conjecture. The universality argument is in part tenuously supported by some 
systematic numerical studies in two and three |^ dimensions, although there is no solid evidence for it 
yet in four dimensions |l^ . The authors of Ref. ||l^ have also independently emphasized the importance 
of exploring the effects of different choices for the measure parameter a. In the context of dynamical 
triangulations and their relationship with simplicial gravity, the irrelevance of the measure parameter has 
also been proposed in . Finally let us mention that one can argue backwards [pj| that the Regge lattice 



measure of Eqs. (6.14) and ( |6.16| ) provides further support for the correctness of the continuum DeWitt 



functional measure approach. 

The above measure can also be obtained by considering a simplicial analog of the DeWitt supermetric. 



as was suggested in M, E8|. One writes for the induced metric of Eq. (2.7) 



6g{s)r = E G''''(.g(s)) S9.,,{s)5gM{s) , (6.17) 



°The functional measure over edge lengths in Eq. ( |6.16[ ) does not have compact support, and a cosmological term (with 
coefficient A > 0) is therefore essential in obtaining convergence of the functional integral for large edge lengths. 



with the inverse of the lattice DeWitt supermetric given now by 



G'^^^'W)] = h^W) W\^)9'\s)+g^\s)g^\s) + \g^={s)g^\s)\ , (6.18) 

and (A 7^ —2/d). This defines a metric on the tangent space of positive real symmetric matrices g. The 
resulting functional measure is the one of Eq. ( |6.1E| ), with, by construction, = (d - 4) (d + 1) /4, the DeWitt 



value. Thus 

j dti[f] = jY{[deiG{9{s))\^\{ d9^J{s) , (6.19) 

s i>j 

with the determinant of the super-metric G^^^^ {g{s)) given by 

detG{g{s)) (X {l + \d\) [g{s)f-^^^''+'^^^ , (6.20) 
and therefore up to irrelevant constants 

[df^if] - rn[^(^)riK ' (6-21) 

with a = (d — 4)((i+ l)/4. The measure factor can be exponentiated and written equivalently as an effective 
action contribution 

" " " (6.22) 



'^^[^'l ^ / n< exp|a^logF(s)| 



Its effect is to suppress or enhance, depending on the sign of a, contributions from small volumes. As such, 
it acts much like a cosmological constant contribution 

expj- A ■ (6.23) 

6.3 Lund Regge Approach 



The previous approach to the functional measure is based on a direct discretization of the continuum 
measure, and leads to a unique local measure over the squared edge lengths (modulo the volume factors), 
in close analogy to the continuum expression. Alternatively, one can try to find a discrete form for the 
supermetric, and then evaluate the resulting determinant. 

In a paper Lund and Regge offered a slightly different approach to the measure problem p5| , in connection 
with the 3 + 1 formulation of simplicial gravity. The idea, recently re-analyzed by the authors of Ref . | |4^ , 
was to obtain a lattice analog of the DeWitt supermetric, by considering the quantity 

\\SlY = E G.Al')Slf6l^ , (6.24) 



where Gij{P) plays a role analogous to the DeWitt supermetric, but now on the space of squared edge 
lengths. One way of constructing the explicit form for Gij{P) is to write the squared volume of a given 
simplex in terms of the induced metric within the same simplex, 



Then compare the expansion of the determinant of the metric in the continuum, 

det (g^J + Sg,j ) = exp Tr log(5,j + Sg,j ) 

= det(5„ ) [ 1 + g'^Sg^j + g'^'Sg.jSgki - g'^'Sg^kSgu + ■ ■ 



to the analogous expansion for the square of the volume of a simplex 



(6.25) 



Slf Sli 



(6.26) 



(6.27) 



Identifying terms of order (^Zf )" in Eq. ( 6.27 ) with terms of order {SgijY^ in Eq. ( 6.26 ) one obtains to linear 
order 

1 dVHP) 



Slf = ^ Jdet{g,,)g^^ 5g, 



and to quadratic order 



W E^SS^^^'^^I = i.^^^)[9''9'''Sg,,5gu- g''g''5g,k5gu 



(6.28) 



(6.29) 



Remarkably, the right hand side of this equation contains precisely the expression appearing in the continuum 



supermetric of Eq. (6.2), for the specific choice A = — 2. After summing over simplices one obtains 



i ^y^dct(5.,(s)) [g'Hs)9'\s)+9'\s)9"'is)-'29''{s)9''\s)]dg,,{s)Sgki{s) ^ Y^G,,{l') 6lf 61^ , 

s ij 

(6.30) 

with 
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(6.31) 



which now determines the matrix Gij{P) appearing in the Lund-Regge metric for deformations in the space 
of squared edge lengths, Eq. ( |6.24 ). The analogy with the continuum expression is brought out more clearly 
when one factors out the volume element, and writes 



- E ns) 



v^{s) ^ dif di^ 
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(6.32) 



The volume factor ambiguity present in the continuum measure is not removed though. As in the continuum, 
different measures on the edge lengths are obtained, depending on whether the local volume factor V{s) is 
included in the supermetric or not. In parallel with Eq. ( |6.9| ) one writes therefore 



The metric in edge length space is again left unchanged by this rewriting, but the measure obtained from 
det G depends on a parameter lu' , and by the usual arguments one obtains a parameterized functional 
gravitational measure 

J d^,[f] = / n VdetG|7')(/2)d/2 (6.34) 

with 

Thus the matrix G^" should be thought of as defining a one-parameter family of measures. 



A somewhat disturbing feature of the Lund-Regge lattice measure of Eq. ( 3.34 ) is that it does not give 



the correct result already in one dimension [p8[ . The one-dimensional action for pure gravity of Eq. (3.46) 



proportional to the length of the curve, is invariant under the local gauge transformations of Eq. (3.47), 



Sin = Xn+l - Xn , (6.36) 

where the Xn's represent continuous parameters defined on the lattice vertices. Any local variations of the 
edges which have the above form (and, we should add, also do not violate the constraint In > 0) obviously 
leave the physical length of the curve unchanged. The invariant measure in one dimension is therefore 

N 



J ,;_n 



which falls precisely in the class of measures encompassed by Eq. (3.15), with a = —1. On the other hand 



the arguments leading to the measure of Eq. (6.34) give P — g to zcroth order, — Sg/g to first order, 

and = to second order, and therefore, since G„m = 0, H^^'^IP = 0! Incidentally, in one dimensions a 
physically motivated invariant distance between manifolds is = [L{1) — L'(Z')]^, which is nonlocal. 

Another somewhat undesirable feature of the Lund-Regge metric is that in general it is non-local, in 
spite of the fact that the original continuum measure of Eq. ( 3.10 ) is completely local. ^ After all, metric 



perturbations and squared edge lengths are linearly related to each other, and locality for one measure should 
translate into locality for the other measure. Non-local contribution to the original (non gauge-fixed) measure 

^After imposing a gauge condition, such as the conformal gauge, the measure can become non-local as in ordinary gauge 
theories. But such a gauge fixing term is only necessary in perturbation theory to remove the zero modes of the action, discussed 
in Sections 3. and 4. Nonperturbatively one would expect that no gauge fixing is necessary, as the effects of the gauge zero 
modes are expected to cancel out in averages of physical quantities, as in ordinary lattice gauge theories [pj. 



seem as unattractive as non-local action contributions. A non-local measure makes it virtually impossible to 
study the theory non-perturbatively. On the other hand it is clear that, for some special choices of uj' and d, 
one does recover a local measure. Thus in two dimensions for lo' = —\ one obtains again the simple result, 

/•oo 

rf^i?'] = / n ' (6.38) 

(which incidentally is non-singular at small edge lengths and represents therefore in this respect an acceptable 
measure). This measure appears therefore, on the basis of purely theoretical arguments, to be as good as 
any other measure with a different lo' . Given the possibility of a choice for w', it would seem natural that 
one should chose its value in such a way that the measure has the simplest form, here the local form without 
any volume factors: after all the physical theory should not depend on the bare parameter u' . 

We should also remark that the appearance of non-local measure contributions (and in particular for 
some values of the parameter w', but not for others) makes one question the initial justification for starting 



in the first place with an expression, such as the one in Eq. (6.18), which is local. In conclusion it appears 
that the local measure of Eq. ( |6.16| ) provides the simplest theoretically justifiable starting point, if not the 
only possible one. 

7 Gauge Fixing and Lattice Conformal Gauge 



Regge's simplicial quantum gravity does not require gauge fixing unless one intends to perform 

a diagrammatic perturbative expansion on the lattice fl^ . In this respect, the situation is completely 
analogous to ordinary gauge theories, and one expects the volume of the gauge group, the diffeomorphism 
group in the case of gravity, eventually to cancel out in the expression for physical averages, 

_ !d^,[P]0{P) exp{-/[/^]} 
^ ' !d^i[P] cxp{-/[P]} ^ ' 

The lattice diffeomorphism zero modes discussed here (see Eq. ( |3.21| ) and subsequent expressions) do not 
therefore in principle pose a problem in nonperturbative studies of quantized gravity, such as the ones 
presented in ^]. Zero modes are automatically taken into account as the measure explores gauge- 
equivalent choices of metrics. One important distinction with the formal continuum theory is the presence 
of a cutoff in orbit space, due to the enforcement of the triangle inequalities. As a result, the gravitational 
functional measure is highly non-trivial. Such a constraint is not seen to any order in the perturbative weak 
field expansion, it is a genuinely non-perturbative constraint. 

On the other hand in two dimensions the continumn theory can be studied by perturbative methods, 
which are most suitably applied in the conformal gauge |^6| . It is the purpose of this section to elaborate 




2^ . Typical edge length distribution V{1) in four dimensions, for the lattice analog of the DeWitt measure a 
Eq. (6. It )), close to the critical point at Gc (from ref. 



on the connection between the continuum and the lattice theory, both being formulated here in a particular 
gauge. Let us first summarize the results in the continuum. The weak field expansion is used, and one sets 
as usual 

gf,:^{x) = Sf^„ + K h^,y{x) . (7.2) 

As there is no small parameter in two dimension to play with, one assumes k <^ 1, expands, and then sets 
K = 1 at the end. The easiest quantity to compute is the "graviton" vacuum polarization due to one massless 
scalar particle, a one loop diagram here. It is given by 

t^iv{p,q) = \[5t,vP-{jP + q)-Pt,{Vv + qv)-Pu{Vti + (liJ.)] ■ (7.3) 

The calculation of the integral is easily done using dimensional regularization or by the methods of 
p7| . In either case one obtains 

n^i.,a/3(g) = ■^^{q^5f,^-qf,qu)^{q^5ai3-qaqf}) ■ (7.4) 

For a D-component scalar field, the above result is simply multiplied by a factor of D, and the effective 
action, to lowest order in the weak field expansion, is then 



f (Pq 

leff = I hf,^{q)n^,,p^{q) hp^{-q) . (7.5) 



In the conformal gauge, coordinates are chosen which are locally orthogonal, so as to bring the metric into 
the form 

5m.(x) = ,5^.e'^(^) . (7.6) 

Then one has for the scalar curvature 

R{q) = {qp.qu-5puq^)hp^{q) = q^ (p{q) , (7.7) 

and one can therefore re-write the effective action in the form 

leffi^) = / (0 q' ^i-q) = / [idp^r + (A - A,) e^] . (7.8) 

On the lattice one can perform a similar computation, using again perturbation theory [ p8[ . The lattice 
Feynman rules are written down, the integration over the scalar is performed, and an effective action results, 
which can be expanded out in the weak field limit. As the scalar couples invariantly to the gravitational 
degrees of freedom, one would expect that the result should be eventually expressible in terms of invariants. 



Indeed in the continuum one can re- write the effective action of Eq. (7.8) in an invariant form, 



i J cPxd'yR^{x){x\^\y)R^{y) , (7.9) 



where is the continuum covariant Laplacian, = d^^/gg^'^d,^. On the lattice this expression has an 
obvious invariant counterpart pl| , 

5 E ^'J^vl../"' ' (7-10) 



hin0;cs/i,/i.' 



which is obtained from the correspondence between lattice and continuum curvatures derived in [Q. A 
here is the nearest-neighbor covariant lattice Laplacian, as obtained from the discrete scalar action (see 



Eqs. ( |5.9| ) and ( 5.48 )); for a recent discussion of the discretization of this term see also ||5^. It introduces 
an effective long-range interaction between deficit angles. For this reason it is actually preferable to study 
non-perturbative aspects of the model leaving the scalar fields un-integrated, which keeps the action local 

The previous discussion provides a background for motivating the introduction of the conformal gauge 
on the lattice. It is legitimate to ask therefore what is the simplicial lattice analog of the gauge condition 



of Eq. (L6). The conformal gauge implies a local choice of orthogonal coordinates. It is clear from the 
discussion in Section 3.1 that there is a corresponding choice on the lattice. Indeed in the development of 
the weak field expansion a uniform orthogonal set of coordinates was chosen, with a diagonal background 
metric (with edge length assignments li — 1 for the body principals {i — 1,2) and l^ — \/2 for the diagonal 
{i = 3)). For these coordinates (see Eq. ( |3.16 )) one has for the background metric 

al? = % • (7.11) 

The lattice conformal gauge choice corresponds to an assignment of edge lengths such that locally 



y ^(liin) - liin) - liin)) li{n) , 

Here the lattice fields ip{n) have to be defined on the lattice vertices, and so are the gauge degrees of freedom 



Xfiin), as can be inferred from Eq. ( 3.23 ). It is clear therefore that a choice of lattice conformal gauge 
corresponds to a re-assignment of edge lengths about each vertex, in such a way that the local curvature is 
left unchanged, but at the same time the induced metric is brought into diagonal form; it corresponds to a 
choice of approximately right-angle triangles at each lattice vertex. 

This result is further illustrated in Figures 22 to 24. The surface shown in Figure 22 has been brought 
into the lattice analog of the conformal gauge, by re-assigning edge lengths in such a way that individual 
triangles look as close as possible to right-angle triangles. In going from Figure 24 to Figure 23, repeated 
gauge transformations must be performed on the vertices, by reassigning edge lengths in such a way that 
local areas and volumes are kept unchanged. It is easy to see that such a construction can always be done, 
except in some rather pathological cases. 



^''it is encouraging that the conformal mode stays massless in the full non-perturbative treatment of the two-dimensional 
simplicial lattice theory, without the necessity of any sort of fine-tuning of bare parameters [ po| , pl| . 



Fig. 21. Description of a smooth surface in the lattice analog of the conformal gauge. 




Fig. 22. Enlarged view of a small region on the surface in Fig. 21. 




Fig. 23. Gauge equivalent description of the enlarged view, of a small region of the original surface, represented in 
Fig. 21. 

The gravitational contribution to the effective action in the lattice conformal gauge can, at least in 
principle, be computed in a similar way. Let us sketch here how the analogous lattice calculation would 



proceed; a more detailed discussion will be presented elsewhere. In the continuum the metric perturbations 
are naturally decomposed into orthogonal conformal and difFcomorphism parts, 

Sg^^ix) = 9p,,{x) S(p{x) + ^pXi^{x) + V^x^(a;) . (7.13) 

where V,y denotes the covariant derivative. It should be clear from the discussion in Section 3. that a 
rather similar decomposition can be done for the lattice degrees of freedom, by separating out the lattice 
gauge transformations (which act on the vertices and change the edge lengths without changing the local 



volumes and curvatures) from the conformal transformations (which do change them) in Eq. ( 3.39 ). The 



explicit form for the lattice diffeomorphisms, to lowest order in the lattice weak field expansion, is given in 



Eq. (3.21), while the explicit form for the lattice conformal transformations in given in Eq. (3.27), which 
makes it obvious that such a decomposition can indeed be performed on the lattice. In the continuum, after 
rewriting the gravitational functional measure in terms of conformal and diffcomorphism degrees of freedom, 

J dfi[g] = J dti[ip]d^i[x][dct{L+L)]^ , (7.14) 

one has to compute Jacobian of the operator L. It is determined, in the continuum, from 

(L+Lx)^ = ^''i^^.X. + ^.X^~9^..^'Xp) , (7.15) 
One then obtains for the effective action contribution in the conformal gauge, 

[det{L+L)y^ ^ exp{-/e//((^)} . (7.16) 



with leff of the form in Eq. (7.5) to lowest order in the weak field expansion. A diagrammatic calculation. 



similar to the one for the scalar field contribution, gives in the continuum the celebrated result g^, |3^ 

13 1 

np^pa(g) = ■^{qpiqu~5puq^)-^{qpq„~5p„q^) . (7.17) 



On the lattice the functional integration is performed over the squared edge lengths, as in Eq. ( 6.15| ) 



But it seems a technically challenging task to compute that Jacobian that maps the edge length variables 



(which define the lattice metric) to the orthogonal lattice diffcomorphism variables of Eq. (3.23) and the 



lattice conformal fields of Eq. ( ^.27 ), with the appropriate Jacobian included. It is also quite possible that 



one might have to go beyond the lowest order in the lattice perturbative expansion. 
As a consequence the total Liouville action for the Liouville field (p — -^R becomes 



967r 

To lowest order in the weak field expansion the critical value of D for which the action vanishes is Dc — 26, 
but this number is modified by higher order quantum corrections. In any case, for sufficiently large D one 



expects an instability to develop. Numerical nonperturbative studies of two-dimensional gravity suggest that 
in the lattice theory the correction is large, and one finds that the threshold of instability moves to values as 
low as Dc ~ 13 It is unclear if this critical value can be regarded as truly universal, and independent 

for example on the detailed choice of gravitational measure (we are referring here to the choice of parameter 
a). 

8 Conclusions 

We have shown in this paper that Regge's formulation of simplicial gravity is endowed with a remnant 
of the continuous local gauge invariance of the original continuum theory. The appearance of zero modes 
corresponding to the diffeomorphisms in the continuum is particularly transparent in the weak field expan- 
sion. Nevertheless, the presence of a local invariance in the discrete action can be exhibited, via the detailed 
explicit calculations presented in this paper, for almost any conceivable choice of background lattice. It was 
shows in particular that the structure of the zero modes corresponds precisely to the discretized form of 
the diffeomorphism transformation law in the continuum. We have underscored the fact that the squared 
edge lengths correspond to the metric components in the continuum, and that such a result is therefore 
hardly surprising (indeed it has been known for some time in some special cases). Explicit calculations also 
show that the gauge and conformal modes are consistently defined as acting on the vertices of the lattice. 
Although our derivations have mainly been restricted to the two dimensional case, where they are more 
transparent and one does not run the risk of drowning is a sea of indices, wa have argued that they have 
general applicability, and in a number of cases we have indicated the structure of the general result. 

Our results have a bearing on the issue of the gravitational measure in simplicial gravity. As the metric 
degrees of freedom in the continuum correspond to the squared edge lengths in the lattice theory, it is clear 
that functional integration in the latter should be performed again over the edge lengths squared. We have 
provided a number of arguments in support of this statement, based on DeWitt's approach to the functional 
measure in the continuum. We have argued that the lattice measure is essentially no less unique than the 
original continuum (DeWitt) measure, with ambiguities restricted to local volume factors, and which most 
likely are not relevant in four dimensions. It is unlikely that further insight into this issue can come from 
analytical work, and it is hoped that future numerical simulations will support this conclusion, for which 
there is already some partial and incomplete evidence. At the end of the paper we have considered the 
introduction of gauge fixing terms in the lattice action, which are needed in order to remove the gauge zero 
modes of the gravitational action in perturbative calculations. Again the situation is similar to what happens 
in the continuum when one performs perturbation theory, where one first separates out the infinite gauge 



Table 1: Critical exponents of random and non-random Ising models. 
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Onsager solution on regular flat lattice 


1.75 


0.125 








1 


Ising spins coupled to gravity!^ H 


1.73(2) 


0.124(3) 


-0.06(11) 




0.98(1) 


Matrix model and CFT 1^ 


1.333... 


0.333... 


-0.666... 


-1.0 


1.5 


Random Ising spins in flat space 10 


1.32(3) 


0.31(4) 


-0.65(4) 


-0.98(4) 


1.46(8) 



volume contribution. As a specific example, we have discussed how one goes about constructing the lattice 
conformal gauge. 

A more practical motivation for our work has been to try to understand the recently discovered discrep- 
ancy between the critical exponents for matter coupled to gravity in two dimensions as computed in the 
lattice regularized model for gravity pO| , |3^ , and the corresponding conformal field theory predictions for 
central charge c = i Particularly significant in this respect appears to be the recent realization 

that the conformal field theory exponents describe two-dimensional random systems in flat space, and do 
not correspond to "gravitational" dressing of correlators |Q (see Table I) . Recent independent calculations 
have to some extent confirmed this result []6l| . 

Previously the authors of Ref. |6^ had considered the Dirac equation on a two-dimensional lattice where 
sites have been removed randomly - a doped lattice. They argued that in this case the fermions acquire a 
quartic interaction and become Thirring fermions, thus changing the critical exponents and the universality 
class. In our opinion, the results for simplicial gravity found in pO| , and the conformal field theory exponents 
of Issj] , can simply be made consistent with each other, if an additive gravitational dressing of critical 
exponents is taken to be zero in both cases, for both non-random and random matter (which are known to 
have different critical exponents in flat space to begin with). 
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